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THE DETERMINATION OF COPPER BY MEANS OF ATOMIC 
ABSORPTION SPECTROSCOPY* 


by 
A. STRASHEIM, F. W. E. STRELOW and L. R. P. BUTLER. 


OPSOMMING 


Die bepaling van koper in ertse, konsentrate en in ertsafal met behulp van atoomabsorpsie, 
word beskryf. ’n Vinnige ekstraksieprosedure om die kopererts ens. in oplossing te bring, met 
behulp van gemengde sure, maak die metode geskik vir produksiekontrole. Vir koperkonsentrasies 
wat varieer van 0-18 tot 24% is die noukeurigheid van die metode beter as 7%. 


SUMMARY 


The determination of copper in ores, concentrates and tailings by means of atomic absorption, 
is described. A rapid extraction procedure to bring the copper ores etc. into solution, using acid 
mixtures, allows the method to be used for production control. The accuracy of the method is 
better than 7% for copper concentrations varying from 0-18 to 24%. 


INTRODUCTION 


The use of atomic absorption spectroscopy for analytical purposes was revived 
by Walsh?! and furthered by Russell, Shelton and Walsh.? Other authors*. * have 
since been equally successful in developing analytical methods based on this principle. 
Atomic absorption spectroscopy, when used for analytical purposes, has certain 
advantages over other methods, viz. : 


(a) It can be used for absolute determinations. } 


(b) It does not require standards in the same physical condition, but does require 
that samples and standards should be in solution. 


(c) It is not affected by self reversal or self absorption. 
(d) The very stable hollow cathode excitation gives good accuracy. 


These advantages suggested its use as a routine method for the analysis of copper 
in sulphide or silicate ores. Methods have been described for the estimation of copper 
by emission flame-photometry,® but they have the inherent difficulties experienced 
with this type of analysis. In this paper a method is discussed for using atomic 
absorption spectroscopy for determining copper in solutions prepared from ores, con- 
centrates or tailings. 


EXPERIMENTAL 


Apparatus. A Hilger and Watts Uvispek spectrometer together with atomic 
absorption attachments, was used. 


The jet in the water-cooled burner was replaced by one having a 0-027 inch 
diameter hole to enable it to be used with Handigas fuel. It was found that if tap 
water was used to cool the burner, this became too cool and vapour condensed on 
the lower face, causing contamination and less efficient vaporisation. These difficulties 
were overcome by circulating water from a 5 litre reservoir through the burner top 


*A paper read in Stellenbosch at the Fourteenth Annual Convention of the S.A. Chemical 
Institute on February 11, 1960. 
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by means of a Stuart-Turner centrifuge pump. During use, the temperature of the 
water rose to approximately 60°C and then remained constant. The burner was 
positioned approximately ? inch below the optical axis. This position was found by 
adjustment of the burner height to obtain maximum absorption. — 


The performance of the spectrometer was appreciably improved by stabilising 
the mains supply to the hollow cathode and Uvispek power supplies with a Wandel 
and Golterman voltage stabiliser. 


Principle of the method. The principle of the method has been fully described by 
Walsh and others }. », 3.4 and will only be briefly mentioned here. The light from the 
cathode passes through the flame into the spectrometer which is set at the wave- 
length of one of the resonant lines of copper, viz. 3247A. If now a solution containing 
copper is sprayed into the flame, sufficient heat is available to break up the solution 
into an atomic vapour. A large number of free neutral copper atoms will then be 
present in and around the flame mantle. These atoms absorb the characteristic 
radiation of copper coming from the hollow cathode lamp, and the spectrometer 
measures the decrease in light intensity at this resonant wavelength. 


General analytical procedure. In the atomic absorption method using the Uvispek 
spectrometer the analytical results are based on single absorption measurements. 
Constancy of the apparatus is, therefore, a stringent prerequisite. New analysis 
curves must thus be prepared when the instrument has to be used for a new analysis 
or a new series of analyses and standard solutions have also regularly to be sprayed 
into the flame to check the instrument. Readily available material for standards, 
which can easily be prepared, is thus a necessity. A primary standard solution for the 
preparation of the analytical curves and for checking the instrument was prepared by 
dissolving a weighed quantity of specpure copper in warm dilute hydrochloric acid. 
This solution was made up to volume to give a standard stock solution containing 
1000 p.p.m. of copper. Further standards were prepared by dilution from this stock 
solution. 


The general analytical procedure was as follows: After allowing sufficient time 
for the instrument to warm up, a series of standards was sprayed through the flame, 
the different density readings recorded and the analytical curve prepared (a typical 
curve is illustrated in Fig. 1). The solutions containing the unknown amounts of 
copper were then sprayed into the flame, the density readings recorded and the 
unknown concentration values read from the analytical curve. 


It was found unnecessary to rinse the spray chamber after each solution, as 
fresh solution removed contamination within a few seconds. 

The operating data of the instrument are given in Table I. The settings were 
found by adjusting all components to optimum conditions. 


TABLE I 

Standard analytical conditions 
Element line | Cu 3247A 
Slit | 0-20 mm 
Meter sensitivity setting 
Air pressure 20 1b/sq inch 
Handigas pressure 0-75 Ib/sq inch 
Sample consumption | 20 ml/min 
Anode current of hollow cathode | 40mA 


(a) Factors affecting atomic absorption results. 


(i) Influence of gas and air pressures. The temperature of the flame is not very 
dependent on the air pressures applied to the atomiser (varying from 15 to 30 Ib per 
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square inch). For air pressures above 30 lb per square inch the angle of the curve is 
lowered as is illustrated for a pressure of 35 Ib per square inch in Fig. 2. Although 
30 lb per square inch gave the steepest curve, a pressure of 20 Ib per square inch was 
chosen as standard operating condition, because of the lower rate of sample consump- 
tion and the more reproducible results. The curves illustrated in Fig. 2 also indicate 
that low concentrations of copper are less affected by changes in pressure. This fact 
should be borne in mind when doing routine analysis. These results also suggest the 
use of two dilutions when determining unknown samples. In this way a check can be 
achieved on the results. 
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Fig. 1—Typical analytical curve 
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In contrast to flame-photometric spectroscopy ® 7 where high temperatures are 
now preferred, the temperature of an atomic absorption unit should only be sufficient 
to vaporise the sample and dissociate the molecules into neutral atoms. When the 
gas pressure to the burner is changed from 0-25 to 1 lb per square inch, the total 
variation possible, the temperature increases from 950 to 1390°C. As no difference in 
the evaporation or dissociation of the copper molecules could be detected with 
variation in gas pressure, 0-75 lb. per square inch was chosen as the standard analytical 
condition. 


(it) Influence of anode current of hollow cathode lamp on analytical curve. Fig. 3 
shows the effect of the anode current of the hollow cathode lamp on the analytical 
curve. The lower the anode current, the steeper and the more linear the analytical 
curve becomes. A decrease in the anode current results in a decrease in the cathode 
temperature. Under these conditions the copper line is narrowed and more effective 
absorption of the copper resonant radiation can be expected, as absorption occurs 
in a very narrow band of wavelengths. 


The line intensity of the hollow cathode lamp is unfortunately related to the 
anode current and the sensitivity of the spectrometer limits the use of anode currents 
which are too low. These findings are in agreement with those of Allan. 


(iii) Influence of extraneous elements and acids. One of the drawbacks of flame- 
spectroscopy is that chemical combinations occur in the solution and in the flame. 
This prevents the atoms from coming into the flame in a free neutral form. In order 
to test the influence of extraneous elements, a series of solutions was prepared con- 
taining varying concentrations of Ca, Na, Fe, K, Al and Co, the elements most likely 
to occur as impurities in copper ores and associated materials. The influence of acids 
was also tested and the results in Table II show that no influence could be detected. 


(b) Procedure for the analysis of copper ores. The success of an analytical method 
when used for production control generally depends on the “lapse of time’’ before 
results become available. A rapid extraction procedure to bring the copper ores etc. 
into solution was thus sought. Various methods of dissolving the samples were tried. 
The most successful acid mixtures to dissolve a 0-5 gram sample were, in order of 
merit: 

(i) 10 ml 40% HF + 10 ml conc. HNO,, 


(ii) 10 ml conc. H,SO, + 10 ml conc. HNO, and 
(iii) 10 ml conc. HCl + 10 ml conc. HNO, 


The general extraction procedure is as follows: The samples are weighed and 
transferred to platinum crucibles. After adding the HF-HNO, mixture, the crucibles 
are heated on a hot plate. The solution is then transferred to a 200 ml volumetric 
flask and made up to volume. This solution is filtered, the first 20-30 ml discarded and 
the next 20-30 ml used for the analysis. If the copper content of this solution is too 
high, the necessary dilutions are made. 


With some flotation samples a little black carbonaceous or sulphurous matter 
remained undissolved after the acid treatment. The amount of copper present in this 
residue was found to be negligible. 


To determine the optimum extraction time, three samples containing 19-6, 
27°5 and 30-64% copper were digested for varying times. The copper concentrations 
in the extracted solutions were then determined using the atomic absorption pro- 
cedure. These results are shown in Table III. From these results it is evident that: 


(i) for full extraction a digestion time of about 40-60 minutes is required, 
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(ii) for routine use of the method an extraction time of 20 minutes may suffice, 
if relative standards are used for the preparation of the analysis curves. 


In the latter case, normal routine samples of known copper content can be used 
as standards. The advantage of this method is that the standards and unknown 
samples are subjected to the same chemical treatment. 


TABLE III 
Influence of extraction time on analytical values. 


Concentration of Cu 


| 
based on mine Extraction time Atomic ab- | % Difference between 
| chemical valuesfor | for atomic ab- (sorption values,| chemical and atomic 
Samples ores, % Cu | sorption, (min) % Cu | absorption values 
426 | 19-6 5 15-5 —21:1 
3677 27-5 13 —52-9 
3685 30-64 22:8 —25-5 
Average —32-2 
| 
426 19-6 10 |} 175 —10-9 
3677 | 27°5 25 — 91 
3685 30-64 | | 27-44 | —10-5 
| | Average —10-2 
426 19-6 | 20 | 18-5 | — 5-7 
3677 27:5 26-1 — 5-2 
3685 | 30-64 28-54 — 6-7 
| | Average — 5-9 
426 19-6 | 40 19-4 — 1-02 
3677 27-5 25-7 — 671 
3685 30-64 : 29-14 — 4-88 
Average — 4-20 
426 19-6 Treatment A 19-6 | — 0-135 
3677 27-5 60 26-6 | — 3-27 
3685 30-64 29-5 — 3-83 
| 
| Average — 2-41 
426 19-6 Treatment B 18-3 — 6-39 
3677 27-5 60 27-9 | + 1-35 
3685 30-64 30-38 — 0-84 
| | Average — 1-9 


Treatment A: Ores dissolved directly and treated for 60 min. 
Treatment B: Ores ignited first in porcelain crucible then treated as in A. 


RESULTS 
From the analytical curve illustrated in Fig. 1, it is evident that the linear 
portion of the curve extends to approximately 30 p.p.m. of copper. Above this value 
accurate results are still obtainable, but as concentrations become higher, they 


become less reproducible. The lower limits of analysis are limited by statistical fluctua- 


tions inherent in the apparatus. 
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From Table IV it can be seen that a concentration of 0-1 p.p.m. of copper can be 
detected. At this concentration the coefficient of variation is, however, above 10%. 
From 1 to 50 p.p.m. the average coefficient of variation is 1-76%, while at about 
100 p.p.m. this value increases to about 3%. 


TABLE IV 
The change of coefficient of variation with different concentrations of copper solution 


Concentration Coefficient of 
of Cu in solution,| variation on 5 
p-p.m. repeat analyses 
0-106 11-8 % 
0-128 14-3 % 
0-182 10-1 % 
0-32 84% 
0-44 97% 
0-68 % 
1-19 3-80% 
1-88 2-31% 
2-4 1:16% 
2-6 2-60% 
3-3 1-20% 
45 1-69% 
9-5 0-80% 
19-03 1-40% 
38-0 1-07% 
47°8 
96-94 291% 


When analysing solutions with high copper concentrations, dilutions have to be 
made. Personal error in measurement is magnified under these circumstances. Very 
concentrated solutions should also be avoided as these gave consistently low readings. 
This is probably due to poor spraying in the atomiser. ? 


The day to day reproducibility of the analytical curve is good, provided that all 
components are adjusted to optimum conditions. 


Information in connection with the accuracy of the method may be gained from 
the results tabulated in Tables V, VI and VII. The chemical values are based on 
results of a standard chemical method ® using the iodide-thiosulphate titration with 
starch as indicator. From these results the following may be concluded: 


(i) Using a 20-minute extraction procedure and relative standards, copper 
concentrations in ores can be determined with an accuracy of 5-3%. This accuracy 
seems satisfactory when compared with the results obtained by the two chemical 
methods. 


(ii) Using absolute standards, low values are recorded for high copper concentra- 
tions. The X-ray fluorescent method, however, gives results which are in agreement 
with the atomic absorption results. The accuracy of the method under these condi- 
tions is +6-2% for copper concentrations varying from 0-18% to 24% copper. 


(iii) The atomic absorption method seems especially suitable for the determina- 
tion of copper in ores with concentrations varying from 0-1% to 5% when absolute 
standards are used. With proper extraction techniques, the results given in Table VII 
seem to indicate that with the use of relative standards, copper concentrations as 
high as 25% can very satisfactorily be determined. 


(iv) The chemical method as used at the mines seems unreliable for low copper 
concentrations. 
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TABLE V 


Comparison of chemical analysis of ores and atomic absorption analysis of solutions obtained with 
20 minute extraction using relative standards. Results given in percentage copper 


| 


| 


C.S.LR. Atomic Percentage difference 
Sample | chemical | absorption | between chemical and 
No. values | values atomic absorption 
(rel. stds) values 
1 24-18 24-18 0 
2 19-58 19-58 0 
3 9-58 8-81 — 8-04 
4 4-16 4-42 + 6-25 
5 4-12 4-00 + 2-91 
6 3-34 3-67 + 9-88 
7 2-70 2-80 -+ 3-70 
8 2-38 2-30 — 3-36 
9 1:95 | 1-95 0 
10 0-81 | 0-76 — 617 
11 | (054 (0595 410-19 
12 0-37 | 0-35 — 5-41 
13 0-26 | 0-285 + 9-62 
14 0-21 | 0-20 — 4-76 
15 0-18 0-17 — 5-56 


TABLE VI 


Comparison of results obtained with different analytical methods. All atomic absorption results were 
obtained with the use of absolute standards. Results are given in percentage copper. Values marked * 
were not included in the calculation of standard deviation 
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| 
| Percentage dif- 
X-Ray ference of 40 Percentage 
| Atomic fluores- } min. extraction difference 
| absorption cence | atomic absorp- between 
Extrac- | values values C.S.LR. Mines tion values from | C.S.I.R. and 
tion } chemical | chemical | C.S.I.R. chemi- | Mines chemi- 
time 20 min | 40 min | 40 min values values cal values cal values 
Sample | | | | 
No. | } 
1 22-6 23-78 22-81 24-18 24:69 | — 1-65 — 2-11 
18-1 18-43 18-0 19-58 19:77 | 5-87 — 0-97 
3 7-63 9-06 9-48 | 958 | 957 | — 5-43 + 0-10 
4 3-59 3-95 35 | 416 415 | — 5-05 + 0-24 
5 3-97 4-01 412 | 3-99 | — 2-67 + 3-16 
6 3-40 3-25 3-34 3-19 — 2-69 + 4-49 
7 2:59 2-57 2:70 26 | — 481 + 3-70 
8 2-12 2-15 238 2:35 | — 9-66 + 1-26 
9 | 1-87 1-82 195 | 1:36 | — 6-67 + 4-62 
10 0-773 0-83 0-81 0-84 | + 2-47 — 3-70 
11 0-564 0-60 0-54 0-69 | +11-11 —27-78* 
12 0-339 0-35 | 0-37 0-45 | — 5-41 —24-32* 
13 0-262 0-25 0-26 | 0-35 | — 3-85 —34-62* 
14 | 0-206 | 0-23 Q-21 0-22 + 9-52 — 4-76 
15 0-177 0-186 | 0-18 0-19 + 3-33 — 5-55 
Percentage standard deviation | 6-23 3-56 


. 
Percentage standard deviation = 5-28 


1960, Vol. XIII JOURNAL OF THE SOUTH AFRICAN CHEMICAL INSTITUTE 81 


TABLE VII 


Comparison of results obtained by chemical analysis with those from atomic absorption using a 
standard chemical extraction method ® 


Atomic 
Chemical | absorption % 
Sample values | values Difference 
1 24-18 23-4 —3-23 
2 19-58 18-6 —5-01 
3 958 | 9-22 | —3-76 
4 416 | 3-99 —4-09 
5 4-12 | 394 


Percentage standard deviation =4-62% 


CONCLUSION 


It is evident from the results given in this paper that atomic absorption spectro- 
scopy provides a relatively quick method for the analysis of copper in solutions and 
in materials containing copper that can be brought into solution. Analysis of copper 
ores, concentrates and tailings can be done satisfactorily for production control 
using a 20-minute extraction.* For accurate analytical results, atomic absorption 
measurements at two different concentration values are recommended. Satisfactory 
results were also obtained using the Zeiss spectrometer with the Hilger Atomic 
Absorption Attachment. 


Thanks are due to Messrs. Hilger and Watts for the loan of the Uvispek spectro- 
meter, to Misses M. C. Maskew and R. E. Rethemeyer, Mrs. C. J. C. Bothma and 
Messrs. E. J. de Villiers and R. P. M. Procter, of this Laboratory, for technical 
assistance; to Dr. A. J. Rossouw, of Rand Mines, Ltd. for his help with the X-ray 
fluorescence results and to Dr. J. N. van Niekerk, of this Laboratory, for the ore 
samples. 
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A STUDY OF THE ANIONIC EXCHANGE BEHAVIOUR OF THE 
VANADIUM(IV) THIOCYANATE COMPLEX. SEPARATION OF 
VANADIUM FROM ALUMINIUM 


by 
L. D. C. BOK and V. C. O. SCHULER 


OPSOMMING 


Vanadium is bevredigend teruggehou toe ’n oplossing van vanadium(IV) in soutsuur, waarby 
tiosianaat-ione gevoeg is, deur ’n kolom van die uitruilingsharpuis Amberlite IRA 400 gestuur 
is. Die adsorpsie en die afwas, asook die geaardheid van die geadsorbeerde vanadiumtiosianaatioon, 
is bestudeer. ’n Analitiese skeiding van vanadium en aluminium is bereik. 


SUMMARY 


Vanadium was satisfactorily retained when a solution of vanadium(IV) in hydrochloric acid, 
to which thiocyanate ions had been added, was passed through a column of the ion exchange 
resin Amberlite IRA 400. Adsorption and elution studies are described, and the nature of the 
adsorbed vanadium thiocyanate ion has been studied. An analytical separation of vanadium and 
aluminium has been achieved. 


INTRODUCTION. 

In the investigation described below an attempt has been made to separate 
vanadium(IV) and aluminium. On adding thiocyanate ions to a solution of vanadium 
(IV) an anionic thiocyanate complex is formed which can be adsorbed on an anionic 
exchange resin, whilst aluminium is not adsorbed iin this way. Teicher and Gordon?! 
have described a method for separating iron(II) ph aluminium by similarly form- 
ing the anionic ferric thiocyanate complex ions and adsorbing these on an anionic 
exchange resin. Vanadium(IV) was chosen for this investigation as vanadium in 
hydrochloric acid solution is usually present in the tetravalent state. ? 


EXPERIMENTAL 
Reagents.—Resin: Rohm and Haas Company—Amberlite IRA 400, effective 
size 0-4-0-6 mm (strongly basic anion exchanger). 
Ammonium thiocyanate: Riedel-de-Haen, A. R. 
Aluminium: B.D.H. Reagent Grade aluminium shavings. 


Vanadium: Hopkin and Williams vanadium chloride (vanadyl) 50% w/w 
(standardised by direct titration with standard potassium permanganate’). 


Hydrochloric acid: A.E. & C.I., chemically pure. 


Exchange columns: The columns used, the treatment of the resin prior to use 
and the method of column operation have been described previously. ¢ 


Analytical methods. Vanadium was determined spectrophotometrically employing 
the method described by Crouthamel e a/.5 Aluminium was estimated by the volu- 
metric EDTA procedure described by Schwarzenbach,® sulphur gravimetrically as 
barium sulphate’ whilst chloride was determined electrometrically. ® 
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RESULTS 

Equilibrium experiments. Preliminary equilibrium experiments were conducted 
to ascertain whether the vanadium(IV) thiocyanate complex could in fact be ad- 
sorbed by Amberlite IRA 400 resin. The procedure adopted was that described by 
Freund and Miner ® who equilibrated 2-0 g of oven dried resin with 100 ml of solution 
carrying the ions under test. The distribution coefficient is defined as:— 

Kd= Amount of metal in resin phase _ ml liquid phase 
~ Amount of metal in liquid phase g dry resin 


In all cases the 100 ml solutions were 0-1M hydrochloric acid containing 1-0 
millimole of vanadyl chloride and the ammonium thiocyanate concentration was 
varied from 0-1M to 2-0M. In addition to aqueous solutions, mixtures of water- 
acetone (1:1) and of water-ethyl alcohol (1:1) were used. The organo-water mixtures 
were investigated because of their application in spectrophotometric analysis where 
the organic solvent increases the optical density of solutions of coloured metal 
thiocyanate complexes, !® presumably by decreasing the dissociation constant of the 
complexes. 


The 100 ml test solutions together with 2-0 g Amberlite IRA 400 resin (chloride 
form) were placed in 250 ml glass stoppered Erlenmeyer flasks and allowed to stand 
at room temperature for three days with intermittent shaking to establish equilibrium 
between the resin and the liquid phases. Suitable aliquot portions were then pipetted 
from the supernatant liquids and the vanadium content estimated. The excess 
thiocyanate was first destroyed by treating the aliquot portions with aqua regia 
and sulphuric acid, followed by evaporation until fumes of sulphuric acid were 
evolved. 


The results obtained (Table I) indicate that the vanadium thiocyanate complex 
is strongly adsorbed from the solvents investigated provided that the thiocyanate 
concentration is kept within the range 0-25M to 1-0M. Furthermore, there appears 
to be no advantage to be gained from aquo-organic solvents. 


TABLE I 


Values for the distribution coefficients of the vanadium(iv) thiocyanate complex in various solvents 
2 g IRA 400 resin, 100 ml 0-1M HCl, 1 millimole VOC, 


Distribution coefficients 


NH,CNS 
molarity | Water Water-acetone | Water-ethanol 
0-10 220 | 340 | 685 
0-25 | 3655 1925 | 2450 
0-50 | 4335 2125 | 2245 
1-00 | 3220 | 1630 | 2245 
2-00 | 2875 | 915 | 2085 


Column experiments. Adsorption of the vanadium(IV) thiocyanate complex on 
30 ml Amberlite IRA 400 resin columns was now attempted using aqueous solutions. 
The effect of the thiocyanate concentration on column adsorption from 0-1M as well 
as from 1-0M hydrochloric acid solutions was studied and the amount of vanadium 
that could be retained by such columns was determined. Columns with an internal 
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diameter of 1-5 cm containing 30 ml of resin, to give column heights of 17-0 cm, were 
employed. 


(a) The effect of thiocyanate concentration. 50 Ml solutions were prepared containing 
10 mg vanadium(IV) and having hydrochloric acid concentrations 0-1M and 1-0M 
respectively and ammonium thiocyanate concentrations 0-5M to 2-0M. These solutions 
were then passed through the resin columns, which were first flushed with 100 ml 
of solution having the same hydrochloric acid and ammonium thiocyanate concentra- 
tion as the respective vanadium test solutions. A flowrate of 4 to 5 ml per minute 
was maintained and the effluents were collected in 250 ml volumetric flasks. The 
columns were washed with successive 20 ml portions of wash solution of the same 
strength as used for the preliminary flushing. Washing was continued until exactly 
250 ml of effluent had been collected. Aliquot portions of the effluents were then 
analysed for vanadium as above. From the results in Table II it appears that vana- 
dium(IV) can be efficiently adsorbed as an anionic thiocyanate complex from hydro- 
chloric acid solution in at least the range 0-1M to 1-0M. 


TABLE II 
Effect of thiocyanate concentration on the adsorption of the vanadium thiocyanate complex 


| 


NH,CNS | V in influent, V in effluent, | Adsorption, 
Mo.arity mg | mg % 
(a) 0-IM HCl 
0-50 | 10-0 0-02 99-8 
1-00 10-0 | 0-03 99-7 
1-50 | 10-0 0-02 99-8 
2-00 | 10-0 0-02 99-8 
(b) 1-0M HCl 
0-50 | 10-0 0-03 99-7 
1-00 | 10-0 | 0-04 99-6 
1-50 | 10-0 | 0-03 99:7 
2-00 | 10-0 | 0-04 99-6 


(b) The amount of vanadium(IV) thiocyanate complex retainable by a 30 ml resin 
column, In order to determine the amount of vanadium thiocyanate complex that 
can be held by a 30 ml column of Amberlite IRA 400 resin, solutions containing various 
amounts of vanadium were passed through such columns. The hydrochloric acid 
concentration chosen was 0-25M, it being felt that lower concentrations would not be 
advantageous when separating vanadium from metals which do not form thiocyanate 
complexes. A thiocyanate concentration of 0-5M was employed as this was found to 
be the concentration at which the vanadium complex is best adsorbed from aqueous 
solution (see Table I). 


By analysing the effluents from the above column experiments (Table III) it 
was found that at least 30 mg of vanadium(IV) can be adsorbed without any signifi- 
cant leakage. The adsorption of larger amounts of vanadium was not attempted 
because the column carrying 30 mg of vanadium had only about 3 cm of clear resin 


(a 
4 
: 
= 
4 
= 


1960, Vol. XIII JOURNAL OF THE SOUTH AFRICAN CHEMICAL INSTITUTE 85 


remaining at the bottom, the upper portion being green in colour due to the ad- 
sorption of the complex. 


TABLE III 


Adsorption of vanadium(iv) thiocyanate complex by 30 ml resin columns from 
0-5M NH,CNS/0-25M HCl 


influent, mg effluent, mg | of 
10-0 | 0-03 | 99-7 
15-0 | 0-04 | 997 
20-0 | 0-03 | 99-9 
30-0 | 0-04 99-9 


' 


(c) Elution of vanadium from the column. Several elutriants were tried in an 
attempt to find a rapid means of removing the adsorbed vanadium complex from the 
columns. The elutriants tested included: water, 0-25M hydrochloric acid, 5M hydro- 
chloric acid, 10M hydrochloric acid, M/10 sodium hydroxide and 2M sodium hydroxide. 


Solutions containing 20 mg of vanadium(IV) in 0-5M ammonium thiocyanate/ 
0-25M hydrochloric acid were prepared and passed through 30 ml resin columns which 
had previously been flushed with 0-5M ammonium thiocyanate in 0-25M hydrochloric 
acid. After the sorption step, the columns were washed with this flushing solution. 
The elutriants were then passed through at a rate of 3 to 4 ml per minute and suc- 
cessive 100 ml portions of the eluates were collected in 100 ml volumetric flasks. 
The vanadium content of the eluates were determined and the results appear in 
Table IV. 


TABLE IV 
Elution of the vanadium(iv) thiocyanate complex with various elutriants 
0-25M 5-0M 10M M/10 2M 
Elutriant: 
H,O HCl HCl HCl NaOH NaOH 
Mg vanadium in eluates 

Ist 100ml 2-24 0-01 0-01 14-8 1-65 11-1 

2nd ,, 3-12 0-25 1-32 3-1 0-68 3-3 

3rd _,, 1-94 1-10 4-40 1-1 0-44 2-9 

4th ,, 0-77 1-72 4-82 0-4 0-41 1-3 

Sth ,, 0-82 2-02 3-34 0-2 0-30 0-8 

6th ,, 0-57 1-58 2-33 0-1 0-35 0-5 

7th ,, 0-39 1-80 1-27 <0-1 0-41 0-1 

8th ,, 0-28 1-39 0-68 <0-1 0-41 <0-1 

Sth ,, 0-20 1-46 0-34 <0-1 0-43 <0-1 
10th ,, 0-11 0-87 0-38 <0-1 0-40 <0-1 
Total 
vanadium 10-44 12-20 18-99 19-7 5-48 20-0 
eluted | 


From the data in Table IV it appears that of the elutriants tested, 10M hydro- 
chloric acid or 2M sodium hydroxide are the most rapid. 2M sodium hydroxide is 
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preferable, as the 10M hydrochloric acid soon decomposes the thiocyanate in the 
eluates with the liberation of sulphur. This decomposition may already occur in the 
resin column and hence retard or even prevent further elution. It is important to 
note that even distilled water will remove some of the vanadium and can thus not 
be used for washing the columns after adsorption of the vanadium complex if this is 
to be quantitatively retained. 


Nature of the adsorbed complex ion. An attempt was made to determine 
the chemical composition of the vanadium(IV) thiocyanate complex ion adsorbed, 
by analysing resin portions loaded with this complex. 50 M1 solutions containing 
50 mg vanadium(IV) in 0-5M ammonium thiocyanate/0-25M hydrochloric acid were 
slowly passed through 0-75 cm internal diameter columns containing 0-5 g resin 
(dry weight, chloride form), which had first been soaked in water and flushed with 
0-5M ammonium thiocyanate/0-25M hydrochloric acid. After the sorption step, the 
resin was washed with 10 ml 0-25M hydrochloric acid and not water as this would 
have removed some of the adsorbed vanadium (see Table IV). The excess wash solu- 
tion was removed by suction. The resin was now quantitatively washed into beakers 
and treated with sodium peroxide, boiled, and acidified with an excess of 5 ml nitric 
acid. The resin was then filtered off and washed with M nitric acid. The filtrates were 
collected and made up to 100 ml. The vanadium, sulphur (converted to thiocyanate), 
and the chloride contents were determined on suitable aliquot portions. Before the 
chloride could be determined, it was found necessary first to convert the vanadium 
present to metavanadate by treatment with sodium peroxide, neutralising with 
nitric acid and finally boiling with sodium bicarbonate. The total chloride capacity 
of the resin was determined according to the method of Kunin!! and found to be 
3-40 milliequivalents per gram of resin (dry, chloride form). The analyses of four 
0-5 g resin portions loaded with the vanadium complex are given in Table V. 


TABLE V 
Analysis of vesin loaded with vanadium(iv) thiocyanate complex 
Mg vanadium/g resin Mg chloride/g resin Mg thiocyanate/g resin 
19-15 9-23 | 223-6 
19-23 9-22 | 223-0 
16-58 10-64 | 215-2 
19-50 9-22 | 232-6 


It has been stated!? that the complex salts of tetravalent vanadium are all 
vanadyl] derivatives and that the chemical co-ordination number of vanadium(IV) 
is six. Hence two anionic vanadium(IV) thiocyanate complex ions are possible, viz. : 
a monovalent ion [VO(CNS),:H,O]— or a divalent complex [VO(CNS),|=. Knowing 
the exchange capacity of the resin used, it is possible from the data in Table V, to 
calculate the number of thiocyanate groups co-ordinated to the vanadium atoms 
for either the monovalent or divalent complex ion adsorbed. Assuming that any of the 
resin exchange capacity that is not taken up by either chloride or vanadium thio- 
cyanate complex ions is taken up by ordinary thiocyanate ions, the number of thio- 
cyanate groups co-ordinated to the vanadium atoms is given by the ratio:— 


moles thiocyanate not adsorbed directly 


moles vanadium adsorbed 


| 
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For stoichiometry this ratio must be an integer. These calculated ratios (Table VI) 
indicate that either of the two complex ions are possible. 


TABLE VI 
Calculated number of co-ordinated thiocyanate groups 


Monovalent complex Divalent complex 


2-87 3-88 
2-86 3-84 
2-86 3-85 
3-26 4-26 
Average 2-96 3-95 


Sidgwick!? states that the complex thiocyanates of vanadyl salts are easily 
formed and are of the type M,{[VO(CNS),]aq. This implies that the water molecules 
are not co-ordinated to the vanadium atom. Furthermore, the stability of co-ordina- 
tion compounds is often affected by the tendency of external ions entering the co- 
ordination sphere, even at the expense of materials already co-ordinated to the metal. 
Thiocyanate is one of the ions quoted !% as having a pronounced tendency to do this. 
It is thus felt that the above considerations favour [VO(CNS),|= for the adsorbed 
vanadium(IV) thiocyanate complex ions. 


Separation of vanadium from aluminium. Vanadium, when present together with 
aluminium, will always partly accompany aluminium while this is being separated 
from other elements. ? It therefore follows that when aluminium is to be determined 
any vanadium present must first be removed. The better known methods to achieve 
this include: 


(1) Precipitation of the aluminium with 8-hydroxyquinoline from an ammoniacal 
solution containing hydrogen peroxide. }4 


(2) Precipitation of the vanadium with cupferron.+® 


The results of the work described above indicated that an anionic exchange 
procedure using a solution containing thiocyanate ions should yield a good method 
for the separation of vanadium from aluminium. 50 M1 solutions containing 10 to 
20 mg vanadium(IV) and 20 to 50 mg aluminium in 0-5M ammonium thiocyanate/ 
0-25M hydrochloric acid were thus prepared and passed through columns of resin as 
described for “‘Effect of thiocyanate concentration”’. The effluents were again collected 
in 250 ml volumetric flasks and the vanadium as well as the aluminium content 
determined. It was found that the aluminium could be determined without first 
destroying the thiocyanate present. The endpoint, however, was more distinct if the 
thiocyanate was first destroyed by evaporation with aqua regia. The latter procedure 
was thus adopted. 


Table VII shows the results obtained when using this method for separating 
vanadium from aluminium. It is clear that good separation of small amounts of 
vanadium, prior to the final estimation of aluminium, is effected. 


rm 
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TABLE VII 
Separation of vanadium(iv) from aluminium 


Influent | Effluent 
Al taken, | V taken, | Al found, V found, 
mg | mg mg | mg 
20-0 | 10-0 | 20-0 | <0-01 
20-0 | 10-0 | 20-0 | <0-01 
50-0 | 10-0 50-2 | <0-01 
50-0 | 20-0 | 50-4 <0-01 
50-0 | 20-0 | 50-2 | 0-02 
50-0 20-0 | 50-2 <0-01 


. Thanks are due to Dr. D. D. Howat, Director, Central Metallurgical Laboratory, 
for making available facilities in this Laboratory, where the major portion of the work 
described above was conducted. 
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THE LIGHT-CATALYSED HYDROLYTIC REACTIONS OF 2,4-DINITRO- 

BENZENESULPHENYL CHLORIDE WITH A NOTE ON THE CARBONYL 

INFRARED ABSORPTION OF ESTER AND AMIDE DERIVATIVES OF 
N-PHENYLHYDROXYLAMINE 


by 
F. KALUZA and G. W. PEROLD 


OPSOMMING 


Dit word getoon dat die hidrolitiese omvorming van 2,4-dinitro-benseensulfenielchloried 
na 2-amino-4-nitrobenseensulfoonsuur deur ’n intramolekulére verskuiwing van suurstof ingelei 
word. Die eersgevormde 2-hidroksielamino-4-nitro-benseensulfoonsuurderivaat word opvolgend 
na die eindproduk reduseer. Albei reaksies kan deur bestraling met lig kataliseer word en vereis 
die inwerking van water. 

SUMMARY 

The hydrolytic transformation of 2,4-dinitro-benzenesulphenyl chloride into 2-amino-4- 
nitrobenzenesulphonic acid is shown to commence with an intramolecular transfer of oxygen to a 
2-hydroxylamino-4-nitrobenzenesulphonic acid derivative, which is then reduced to the final 
product. Both reaction sequences are subject to catalysis by light radiation and involve water as 
reactant. 


2,4-Dinitrobenzenesulpheny] chloride is a very useful reagent for the characteri- 
sation of various types of organic compounds.! The reagent also enters into rapid 
reaction with a variety of small molecules. 


The sulpheny] chloride moiety in 2-nitrobenzenesulpheny] chloride was long ago 
shown? to react rapidly with water to give the sulphenic acid anhydride, while 
aqueous alkali led to the formation of the corresponding disulphide and disulphoxide. 
The intermediate formation of the unstable free sulphenic acid was postulated to 
explain the mechanism of these reactions. Exactly parallel behaviour was subse- 
quently demonstrated for 2,4-dinitrobenzenesulpheny! chloride. * 


A very interesting reaction, involving a larger portion of the reacting molecule, 
was at the time found for both the 2-nitro- and 2,4-dinitro-benzenesulpheny] chlorides 
respectively. Boiling in aqueous methanolic solution converted? 2-nitrobenzene- 
sulphenyl chloride (I) into a mixture of orthanilic acid (III) and 2-nitrobenzene- 
sulphinic acid (IV) in moderate yield. The decomposition of the free unstable 2-nitro- 
benzenesulphenic acid (II) was thought to explain this reaction. 


SCl SOH SO,H SO,H 
NO, : NO, NH, NO, 
I II Ill IV 


Parallel behaviour was later observed in the hydrolytic change of both 2,4-dinitro- 
benzenesulphenyl chloride (V) and the corresponding methy!] sulphenate (VII) to 
2-amino-4-nitrobenzenesulphonic acid (VI) in acid medium, although no sulphinic 
acid was then found. 


c 


a 
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SCl SOCH, 
H*, H,O 4 
(MeOH) (MeOH) 
ON NO, 

V VI VII 
BzCl, 
pyridine 
SO,H 


The yield of the aminobenzenesulphonic acid ((III or VI), respectively) could, 
according to the mechanism proposed by Zincke, * not exceed 50%. The higher yields 
obtained from the related reaction on methyl 2,4-dinitrobenzenesulphenyl chloride 
(V) require a re-interpretation of this transformation. 


Our studies on irradiated solutions of the sulpheny] chloride (V) have now shown 
that this reaction may be catalysed by exposure to light, and that further steps 
are involved, the presence of water still being essential. 


A solution of 2,4-dinitrobenzenesulphenyl chloride (V) in commercial glacial 
acetic acid (containing 0-14°%, water) underwent no change when kept in the dark at 
20-30° for 60 days, the reagent being quantitatively recoverable. When, however, a 
solution of the reagent in acetic acid containing 0-14-2-5% water was exposed to full 
sunlight, or to ultraviolet radiation from a quartz lamp, the colour soon darkened, 
while a crystalline precipitate began to form after a few hours. The product obtained 
after exposure for 70 hours was shown to be 2-amino-4-nitrobenzenesulphonic acid 
(VI), identified by comparison of its S-benzylisothiuronium salt with authentic 
material, by its ultraviolet absorption spectrum, as well as by its reaction with bromine 
water to yield 3-nitro-2,4,6-tribromoaniline* by substitution and by nucleophilic 
replacement of the sulphonic acid group. The 2-amino-4-nitrobenzenesulphonic acid 
was furthermore characterised as its N-benzoate (VIII) pyridinium and S-benzyl- 
isothiuronium salts. 


When a solution of the sulphenyl chloride (V) in anhydrous acetic acid was kept 
in a quartz flask in full sunlight for 35 hours, no precipitate whatsoever was obtained. 


When the acetic acid solution of .2,4-dinitrobenzenesulphenyl chloride was 
irradiated for relatively short periods, a different crystalline product was obtained 
in yields of up to 64%. Although it could not be obtained completely pure, this product 
was identified as the unknown 2-acetoxyamino-4-nitrobenzenesulphonic acid (IX), 
showing strong infrared absorption at 1795 cm-! (C = O in acetoxy).* It was charac- 
terised as the S-benzylisothiuronium salt. 


*A number of reference compounds, derived from N-phenylhydroxylamine, were synthesized 
to allow the decisive identification of ester and amide carbonyl infrared absorption peaks in this 
study. The relevant data are summarised below in the Note appended to this paper. 


O,N NHBz 
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-The new irradiation product (IX) was hydrolysed to acetic acid and the hitherto 
unknown 2-hydroxylamino-4-nitrobenzenesulphonic acid (X).* This compound 
no more showed acetoxy carbonyl absorption in the 1800 cm-! region, but now dis- 
played hydroxyl group absorption at 3279 cm-1!. It was characterised as its dibenzoate 
(XI) pyridinium and S-benzylisothiuronium salts. 


The role of 2-acetoxyamino-4-nitrobenzenesulphonic acid (IX) as an intermediate 
in the change leading to 2-amino-4-nitrobenzenesulphonic acid (VI) was confirmed 
by irradiating a suspension of (IX) in acetic acid in the presence of hydrochloric 
acid when the final product (VI) was readily obtained. 


SOH| SO,H SO,H 
NO, | O,N NHOAc O,N NHOH 
XI IX x 
| 
H,O H,O | BzCl, 
| HCl | pyridine 
| | ACOH 
V 
SO,H SO,H 
AN 
O,N NH, O,N N—OBz 


VI XI 


The overall reaction therefore involves an intramolecular rearrangement, subject 
to catalysis by light radiation and involving water as a reactant, of intermediately 
formed 2,4-dinitrobenzenesulphenic acid (XII). This leads to 2-acetoxyamino-4- 
nitrobenzenesulphonic acid (IX), which is finally reduced photochemically to 2- 
amino-4-nitrobenzenesulphonic acid (VI). This reduction involves water as a re- 
actant, as no smooth reaction is obtained in anhydrous medium. The water presumably 
also accepts the photolysed oxygen atom; while the irradiated solutions do liberate 
free iodine from potassium iodide solution, this does not allow a clear decision as to 
the fate of the photolysed oxygen, because (impure) water is itself known ® to yield 
hydrogen peroxide on irradiation. 

The intramolecular change of (XII) to (IX) recalls the well-known Ciamician 
photochemical transformation® of o-nitrobenzaldehyde (XIII) to o-nitrosobenzoic 
acid (XIV), which proceeds smoothly even in the absence of solvents., 


*Potentiometric titration of this compound in aqueous solution showed pK, = 2-3 + 0-1 and 
pK, = 7:3 + 0-1, clearly reflecting the strong sulphonic acid group and a weak acid group re- 
spectively. As the infrared absorption spectra of compound (IX) and the salts of its N,O-dibenzoyl- 
derivative consistently support the presence of two free adjacent hydrogen atoms on the aromatic 
ring, this suggests that an arylhydroxylamine—>aminophenol change occurs under the conditions 
of this titration. This point is presently being studied. 


| 
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CHO COOH 


XIII XIV 


This similarity was demonstrated by the finding that the reaction studied here 
also occurred on exposing the solid solvent-free reagent in (moist) air to sunlight. 
Many similar examples of light-induced intramolecular transport of oxygen have 
been collated in a recent review.’ These without exception describe transfer® of 
oxygen from a nitro group to an acceptor group situated in an adjacent position on 
an aromatic ring. Furthermore, the acceptor group in nearly all these instances is 
at such a level of oxidation that the reaction stops after transfer of one atom of 
oxygen, with production of o-nitroso compounds. 


In the case of the sulphenyl compounds here studied, however, the sulphur atom 
is at so low a level of combination that it can accept two atoms of oxygen to form a 
sulphonic acid derivative. The nitrogen atom in the original o-nitro* group is in this 
case therefore reduced to a stage beyond that of the o-nitroso group in the instances 
cited. 


EXPERIMENTAL 


M.p.s. are corrected. Infrared absorbtion spectra were taken on the Perkin 
Elmer Model 21 spectrophotometer, in potassium bromide dispersion except 
where shown otherwise. 


2-Amino-4-nitrobenzenesulphonic acid (VI). The sulphenyl chloride (V) * (10-0 g) 
in acetic acid containing 2-5°% water (200 ml) in a quartz flask, was exposed to full 
sunshine for 70 hours at temperatures up to 30°; the precipitate formed was filtered 
off daily. This material was extracted with water to give an insoluble residue (pre- 
sumably bis-2,4-dinitrophenyl disulphide ) (0-29 g) and a water soluble acid (VI) 
(5-39 g, 58% yield). This was recrystallised from water as greyish small flat prisms 
decomposing towards 300° (Found: C, 33-0; H, 2-9; N, 13-0; S, 14-9%; Eq. wt., 221. 
Calc. for C,H,N,O,S: C, 33-0; H, 2-8; N, 12-8; S, 14-7%; Eq. wt., 218); Amax in 
ethanolf 238mzy (log ¢, 4:3), Amin 325 my (log e, 2-9), Amax 390 my (log e, 3-3). 


The acid (VI) (1-95 g) in hot 0-4N NaOH aqueous solution (25 ml) yielded a 
canary yellow sodium salt of low solubility (2-15 g) which was converted?!® to the 
S-benzylisothiuronium salt (3-4 g). Recrystallisation from aqueous ethanol gave the 
thiuronium salt as yellow scales, m.p.{ and mixed m.p. 147-5° (Found: C, 43-8; H, 
4-2; N, 14-6; S, 16-7. Calc. for C,4H,,N,sO;S,: C, 43-7; H, 4:2; N, 14-6; S, 16-7%). 
The reference compound was similarly obtained from the nitroaminosulphonic acid 
(VI) prepared?! from 3-nitroaniline. It melted at 147-5° (Found: C, 43-7; H, 4-2; 
N, 14-6; S, 17-0%). 


*The fact that only the o-nitro group is involved in this reaction supports the intramolecular 
mechanism proposed. 

tCorresponding values for 3-nitro-aniline were Amax 235 my (log e, 4-2), Amin 318 my (log e, 
2-8), Amax375 my (log €, 372). - 

{Kharasch ef al.* report a m.p. of 136-137° for this salt. While in our hands a double m.p. 
was sometimes observed (m.p. 135°, followed by resolidification and m.p. 147°), most of our 
preparations finally melted at 147-5°. 


a\/ hy Y 
| 
ON 
NO, NO 
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The acid (VI) (235 mg) warmed with bromine water‘ gave a precipitate (390 mg, 
m.p. 101-5°). Recrystallisation from glacial acetic acid afforded 3-nitro-2,4,6-tri- 
bromoaniline (yellow needles, m.p. and mixed m.p. 102°) (Found: C, 19-3; H, 0-8; 
Br, 64-1; N. 7-3. Calc. for CgH,Br,N,0,: C, 19-2; H, 0-8; Br, 64-0; N. 75%). The 
reference compound, prepared directly'? from 3-nitroaniline!*, had m.p. 102° 
(Found: C, 19-4; H, 0-7%). 


Salts of 2-benzamido-4-nitrobenzenesulphonic acid (VIII). The acid (VI) (490 mg) 
on treatment with pyridine (2 ml) and benzoyl chloride (2 ml) gave the pyridinium 
salt of the nitrobenzamidosulphonic acid (VIII) (517 mg). Crystallisation from aqueous 
ethanol afforded short yellowish prisms of the pyridinium salt, m.p. 223-224° (Found: 
C, 53-7; H, 3-7; N, 108; S, 8-0. Calc. for C,,H,;N,;0,S: C, 53-9; H, 3-8; N, 10-5; S, 
8-0%); Ymax 1678 cm- 1 (amide C = 


The S-benzylisothiuronium salt was obtained from the pyridinium salt in good 
yield?° and was recrystallised from ethanol-water (1 : 1) as colourless needles, m.p. 
176° (Found: C, 51-5; H, 4-3; N, 11-8. Cale. for C,,H,9N,O,S,: C, 51-6; H, 4:1; N, 
11-5%). 


2-Acetoxyamino-4-nitrobenzenesulphonic acid (IX). The sulphenyl chloride 
(V)(2-0 g) in glacial acetic acid containing 0-14% water* (200 ml) in a quartz flask was 
exposed to full sunshine for 15 hours at temperatures up to 29°. The precipitate was 
filtered off, washed thoroughly with glacial acetic acid and benzene and then dried 
to give the impure acetoxyamino-nitrosulphonic acid (IX) (2-78 g, 59% yield) as 
small scales decomposing towards 300° (Found: C, 33-3; H, 2-8; N, 11-0; S, 125%; 
Eq. wt., 255. Calc. for C,H,N,0,S: C, 34-8; H, 2-9; N, 10-2; S, 11-6%; Eq. wt., 
276); ¥max1795 (C = O in acetoxy), 875 and 832 cm=! (resp. one and two adjacent 
free H on aromatic ring).!*4 


The S-benzylisothiuronium salt was obtained!® by neutralising and reacting 
with S-benzylisothiuronium chloride in ice-cold solution. The product was re- 
crystallised from aqueous ethanol as yellowish crystals, m.p. 123° (Found: C, 43-3; 
H, 4-4; N, 12-7; S, 14-8. Calc. for C,,H,,N,O,S,: C, 43-4; H, 4-1; N, 12-7; S, 145%); 
Vmax 1748 (C = O in acetoxy), 893 and 833 cm-! (1H and 2H as above.) 


2-Hydroxylamino-4-nitrobenzenesulphonic acid (X). The crude acetoxyamino 
compound (IX) (1-0 g) after heating in hot ethanol-water (1:1) (30 ml) gave the 
hydroxylaminosulphonic acid (X)(670 mg), which was recrystallised from water 
yielding small lightly coloured scales decomposing towards 300° (Found: C, 30-7; 
H, 2-6; N, 12-3; S, 13-7. Calc. for CsH,N,O,S: C, 30-8; H, 2-6; N, 12.0; S, 13-7%). 
The infrared absorption spectrum no longer showed ester-carbonyl absorption near 
1800 cm-! but now showed hydroxyl absorption at 3279 cm=}, as well as peaks at 
896 and 884 cm-! (1H) and 852 cm-! (2H). 


The acetic acid obtained in the foregoing hydrolysis was isolated in a separate 
experiment. Crude acetoxyamino compound (IX)(2-28 g) was powdered, dried at 
72°/0-1 mm pressure for 2 hours and then refluxed in 100 ml water until solution was 
complete. The solution was cooled and the crystallised hydroxylamino compound 
(1-89 g) filtered off. The filtrate was neutralised with 0-1 N sodium hydroxide (70 ml), 
evaporated and reacted with p-bromophenacyl bromide (2-0 g) in ethanol (50 ml), 
water (15 ml) and a trace of hydrochloric acid. The product was chromatographed 


*The water content was determined by Karl Fischer titration. 


+The low eq. wt. may be due to the ready hydrolysis (see below), during the titration, of 
the acetoxyamino compound to acetic acid and the hydroxylamino compound. 


¢ 
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over silica gel (100 g) to afford p-bromophenacyl acetate (701 mg, 33%, yield), m.p. 
and mixed m.p. 85° (Found: C, 46-5; H, 3-5. Calc. for C,gHgBrO, : C, 46-7; H, 35%). 


2-N ,O-Dibenzoylhydroxylamino-4-nitrobenzenesulphonic acid (XI) salts. The 
hydrolysed irradiation product (X)(200 mg) in pyridine (1 ml) and benzoyl chloride 
(0-8 ml) gave the N,O-dibenzoate (XI) pyridinium salt (162 mg) which was recrystal- 
lised from ethanol to afford yellowish needles, m.p. 232° (Found: C, 57-3; H, 4-0; N, 
8:2; S, 6-3. Calc. for C,;H,gN,0,S: C, 57-6; H, 3-7; N, 8-1; S, 62%); ?max 892 (weak) 
(1H) and 848 cm-! (2H). 


The pyridinium salt (as above, from 140 mg (X)) was shaken in benzene solution 
with 5% aqueous sodium carbonate solution to yield an insoluble sodium salt (195 mg), 
converted to the S-benzylisothiuronium salt (261 mg) as before. Crystallisation 
from aqueous ethanol afforded the pure salt as yellowish needles, m.p. 212-213° 
(Found: C, 55-5; H, 4-0; N, 9-6; S, 10-8. Calc. for C,,H.4N,O,S,: C, 55-3; H, 4-0; 
N, 9:2; S, 105%); max 844 (2H). 


Conversion of 2-acetoxvamino-4-nitrobenzenesulphonic acid (IX) into 2-amino-4- 
nitrobenzenesulphonic acid (VI). A suspension of the acetoxyamino compound (IX) 
(1-21 g) in glacial acetic acid (50 ml) containing conc. aqueous hydrochloric acid* 

‘ (1 ml) was exposed to full sunshine for 40 hours in a quartz flask. The temperature 
of the suspension rose to 55°. The precipitate was filtered off and dried (1-15 g). It 
(318 mg) dissolved in 3 ml warm aqueous 0-5 N sodium hydroxide solution to yield 
an insoluble red sodium salt (237 mg, 65% yield on IX)) which was converted as 
before into the yellow S-benzylisothiuronium 2-amino-4-nitrobenzenesulphonate 
(above), m.p. and mixed m.p. 146-147° (Found: C, 43-9; H, 4-3. Calc. for C,,H,,.N, O;- 
S,: C, 43-7; H, 4:2%). 


Solvent-free conversion of 2,4-dinitrobenzenesulphenyl chloride (V) into 2-amino-4- 
nitrobenzenesulphonic acid (VI). The dry powdered sulphenyl chloride (V) (1-0 g), 
spread on a watchglass and kept in the sun for 8 hours during one week, turned 
brown and sticky while hydrochloric acid fumes were given off. Extraction with 
benzene left an insoluble portion (480 mg), soluble in water with acid reaction and 
giving a sodium salt (180 mg) which was converted into S-benzylisothiuronium 
2-amino-4-nitrobenzenesulphonate (176 mg), m.p. and mixed m.p. 146-147-5°. 


NOTE 


The carbonyl absorption in the infrared spectra of ester and amide deriva- 
tives of N-phenylhydroxylamine. 


A series of acylated N-phenylhydroxylamine derivatives was prepared, in order 
to check the accuracy of ester and amide carbonyl infrared absorption allocations 
in this field. 


N-Phenylhydroxylamine, m.p. 83° (lit.1® 81°), was converted to the N,O-diacetyl 
derivative, m.p. 44-45° (lit. 1® 43°) which was selectively hydrolysed to the N-acetyl 
derivative, m.p. 68° (lit.17 67-67-5°). 


N-Acetyl-N-phenylhydroxylamine was tosylated in the presence of pyridine to 
yield N-acetyl-O-tosyl-N-phenylhydroxylamine, m.p. 135°. N-Tosyl-N-phenyl- 


~ *Hydrochloric acid was added to simulate the conditions of the overall reaction of (V)—>(V 1). 
The reduction of (IX)—>(VI) also occurred in glacial acetic acid only, though the yield of the 
sodium salt of (VI) was then lower. 
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hydroxylamine, !§ m.p. 143-5° (lit.1® 143-5°) (Found: C, 59-3; H, 5-1; N, 5-1. Cale. 
for CysH,3;NO,S: C, 59-3; H, 5-0; N, 5-3%), was acetylated with acetic anhydride in 
the presence of pyridine to yield N-tosyl-O-acetyl-N-phenylhydroxylamine, m.p. 
138°. N-Tosyl-N-phenylhydroxylamine with benzoyl chloride and pyridine yielded 
N-tosyl-O-benzoyl-N-phenylhydroxylamine, m.p. 155°. 

N-Benzoyl-N-phenylhydroxylamine, m.p. 122° (lit.?° 121-122°) (Found: C, 73-1; 
H, 5-3. Calc. for C,,H,,NO,: C, 73-2; H, 5-2%), was tosylated as above to yield 
N-benzoyl-O-tosyl-N-phenylhydroxylamine, m.p. 111°. 


These compounds showed strong infrared absorption bands in the 1700-1800 
cm-! region; these are tabulated below. 


TABLE 
R, 


Derivatives of N-phenylhydroxylamine, Ph-N-OR, 


% C % H % N 

R,* | R,*| in | in | Mp Formula f 

| ester | amide | | Found | Calc. | Found | Calc. | Found Calc. 

| | | | 
H | Ac — 1645¢ | 68° | C,H,NO, 63-5 63-6 6-0 6-0 9-3 9-3 
Ac Ac | 1799t| 1692t C,9H,,NO, 62-4 62-2 5-8 5-7 7-4 | 7:3 

| | | 
Tos | Ac 7 1695+ 135° | CisHisNO,S | 59-0 59-0 5-0 5-0 4-9 46 
Ac |Tos| 1799 | — | 138° | CisHysNO,S | 58-9 59-0 5-0 5-0 4:7 4-6 
Tos | Bz | — | 1675 | 111° i attaitilidl 65:3 65-4 4-9 4-7 3-9 | 3-8 
Bz |Tos| 1773 155° CyoH,;NO,S| 65-6 65-4 4:8 4-7 3-9 | 3-8 

| 


*Ac = CH;.CO—; Bz = C,H;.CO—; Tos = p-CH,.C,H,.SO,— 


{These infrared spectra were taken in solution in chloroform, all others in KBr dispersion. 
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GAS-LIQUID CHROMATOGRAPHY IN CAPILLARY COLUMNS. 
PART I—THEORY 


by 
P. C. HAARHOFF and VICTOR PRETORIUS 


OPSOMMING 
’n Streng algemene behandeling van liniére nie-ideale gas-vloeistofchromatografie in kapillére 
kolomme word gegee. Alle kolomparameters behalwe die kolomdimensies en die dikte van die 
vloeistoflaag word beskou as kontinue tydsafhanklike funksies van die afstand vanaf die kolomin- 
laat. Daar word aangetoon dat die resultate geldig is mits die kolom ’n groot aantal teoretiese 
plate bevat. 
SUMMARY 


A rigorous general treatment of linear non-ideal gas-liquid chromatography in capillary 
columns is presented. All the column parameters other than the column dimensions and the thick- 
ness of the retentive liquid layer are treated as continuous time-dependent functions of the distance 
from the column inlet. The results are shown to be valid provided only that the column contains 
a large number of theoretical plates. 

INTRODUCTION 

The use of coated capillary columns for gas-liquid chromatography is now 
firmly established.! A study of such columns is important firstly because they appear 
to give higher efficiencies than do packed columns,? and secondly because they are 
more amenable to a rigorous theoretical study than are packed columns. Furthermore, 
a thorough understanding of the processes in coated capillary columns would provide 
a valuable pointer for the development of approximate theories of the other forms of 
chromatography in which the mechanism of the process is more complex. 


The technique may be briefly described as follows. A suitably small gaseous 
sample is injected into an inert carrier gas flowing through a capillary tube on the 
walls of which a thin layer of liquid is deposited. The various components of the sample 
are distributed between the gaseous mobile and liquid stationary phases. Because 
their distribution coefficients will normally differ, the components will be selectively 
retarded by the liquid and they will form bands which move through the column 
at different rates. The efficiency of the separation is impaired by the fact that the 
bands broaden during their passage through the column. It is generally conceded® 4 
that band broadening may be ascribed to two distinct causes—longitudinal diffusion 
of the solute molecules in the mobile and stationary phases, and the resistance of the 
solute molecules to mass transfer in the mobile and stationary phases. As a result 
of the resistance to mass transfer the density of the solute at various points along the 
column cross section will deviate from the equilibrium distribution. This causes the 
solute to move ahead of the band centre when there is an excess of solute in the radial 
region of high flow velocity, while a retarding effect is produced by an excess of 
solute in a low flow velocity region. Longitudinal diffusion of the solute becomes more 


j 
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pronounced as the diffusivity of the solute in the phases increases and also as the 
flow velocity is decreased. The effect of the resistance to mass transfer, on-the other 
hand, is diminished by rapid solute diffusion and slow flow rate. It is evident that the 
highest column efficiency will be obtained when the operating conditions are so chosen 
that both the longitudinal diffusion and the resistance to mass transfer are simul- 
taneously minimised. 


A detailed study of longitudinal diffusion and of the resistance to mass transfer 
in coated capillary columns has been carried out by Golay. 4 An expression is obtained 
in which the column efficiency, in terms of the height equivalent per theoretical 
plate, is related to the various column parameters, e.g. gas flow rate, distribution 
coefficient, column diameter, etc. In the derivation of this question it has been 
assumed that the column parameters are constant at all points along the length of 
the column and also that they do not vary during the time taken by the solute band 
to pass through the column. In practice, these assumptions are generally approxi- 
mately true. Recently, however, a certain amount of attention has been focussed 
on two techniques where the Golay original equation no longer holds. In the first 
technique very long fine bore capillary tubes are employed.? The pressure drop 
across such columns may be appreciable and since the mobile phase is compressible 
the gas velocity may vary significantly along the column length. Several attempts 
have been made to take the variation in the gas velocity into account. 4:5. ® 7,8 
The most elaborate of these has been carried out by Giddings e¢ al. 5. ® who adopted a 
statistical approach. These authors found that their results were not in good agree- 
ment with experimentally determined values.* Golay* has discussed the effect of 
pressure gradients in coated capillary columns in which the resistance to mass transfer 
in the liquid phase is negligible, and the gas velocity is at the optimum value. The 
method used by Golay is roughly equivalent to that of Steward e¢ al.,5.® but it lacks 
a detailed physical background. 


The second technique for which the Golay equation does not apply is one in 
which the column temperature is programmed, ® i.e. the column temperature is pro- 
gressively changed during the passage of the solute bands through the column. 


In the investigation described below an attempt has been made to develop a 
mathematically rigorous theory of linear non-ideal gas-liquid chromatography in 
coated capillary columns. The theory is more general than that given by Golay in 
the sense that all the column parameters, other than the column dimensions and the 
thickness of the stationary phase, are treated as continuous time dependent functions 
of the distance from the column inlet. 


1. Model and basic equations used in present study. Consider a capillary of radius 
T9*, the walls of which are coated with a thin uniform liquid layer of thickness w. 
The radius of the portion of the column occupied by the mobile phase is given by 


* The symbols used are defined on page 113. 
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Under the conditions normally encountered in practice the flow of the carrier 
gas may be assumed to be laminar.!° The carrier gas velocity is then given by 


2 


If a gaseous sample is now instantaneously introduced into the mobile phase 
at the column inlet, the mass flux of the sample in the mobile phase is given by! 


The equation of continuity is 

dc, 

at 


It has been assumed that the sample is sufficiently small so that the distribution 
isotherm will be linear. Equations (2), (3) and (4) may be combined to describe mass 
transfer in the mobile phase: 


1 dc dc dc r?\0 
+ Ox ot s* ax ©) 


The ratio of the diffusivity D of the solute in the gaseous phase to the diffusivity 
D, of the solute in the liquid phase is normally greater than 104. Longitudinal dif- 
fusion in the liquid phase may thus be neglected. 


The mass flux of the solute from the mobile to the stationary phase is given, 
from equation (3), by 


dc 
©) 
and from conservation of mass, by 
dC, 
Ie 


as mass flow in the stationary phase is assumed to occur only in a radial direction. 
The ratio of the volume of the stationary phase to that of the mobile phase is given by 


From equations (6), (7) and (8) it follows that 
2D =) ac, 
Ss r=s 2 ot 


Equation (9) describes the mass transfer of the solute across the gas-liquid interface. 

The following equation, describing the mass transfer of solute into a flat liquid 
layer, has been derived by Golay*: 
-_T ac, 


— Ce (10) 


where 
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and 


An equation similar to equation (10) has been shown to hold for diffusion into spherical 
particles, '? provided that equation (11) is suitably modified to make allowance for 
the curvature of the surface. In coated capillary columns both equations (10) and 
(11) may be used without introducing a significant error provided that w/ry is smaller 
than about 0-3. 


(2) General differential equation describing band development. The final form of the 
required solution of equations (5), (9) and (10) is an equation expressing the average 
concentration of the solute in the mobile phase, ¢,, as a function of time at the column 
outlet. As a first step, equations (5), (9) and (10) are reduced to a differential equation 
containing only one dependent variable, with x and t as independent variables. In 
the general case where the distribution coefficient K varies along the column, it 
becomes mathematically inconvenient to use ¢, as dependent variable, since the total 
amount of the solute in the mobile phase changes during the passage of the band 
through the column. We therefore define a convenient parameter 


where the average concentration of the solute in the mobile diaaat is given by 


From equations (8) and (13) it follows that the total mass of sample per unit length 
of the column is related to c by 


For the purpose of expressing all the dependent variables in terms of c, we begin 
by defining two terms which constitute a measure of the local departure from equili- 
brium: 


8, — ee ee ee ee ee ee ee ee ee (16) 

It is seen from equations (14) and 116) t that 


The ratio of the mass of the solute in the stationary phase to that in the mobile phase, 
is given at equilibrium by 


It follows from equations (13), (17) and (19) that 
= (c — ké,)/(1 + k) an ae. 

and 


5 
P 
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From equations (5), (9), (10), (16) and (17), 


2\ 9 
s*/ 0x 
2D *) acy 23 
a 
ac, 

— (8,),-, = — TK — 24 
at ( ) 

On applying the operator A [ rdr_ to each term of equation (22) to obtain an 


average value, we find by partial oe when use is made of equation (18), 


+2 (om) 
4 r2\ 0 


The boundary flux term may be eliminated by subtracting equation (23) from equa- 
tion (25), and then 7 equation (13) to give 


p? 

It is now convenient to define 


Although equation (27) is introduced for mathematical convenience it will later be 
shown that the value of V at the centre of mass of the solute band is equal to the 
velocity of the centre of mass of the solute band. By substituting equations (20) 
and (27) in equation (26), an equation in terms of c, 8,, and 5, is obtained: 


In order to estimate the relative importance of the various terms, certain approxi- 
mations are introduced. These are later verified by equations (49) and (50), and in_ 
paragraph 6. Following Golay,‘* we shall assume that both 8, and 4, are negligibly 
small in comparison with ¢;, Kc, and c/(1 + k), ice. 


In equations (76), (81) and (83) it is proved that the concentration functions are of a 
Gaussian form; 8,, 52, Cj, 2, c and their various derivatives will thus vary from their 
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maximum to their minimum values in approximately the same distance, at constant 
time, and in the same time, at constant position. It then follows that any derivatives 
of 5, and 6, may be neglected, as a first approximation, in comparison with the 
corresponding derivatives of ¢;, Kc, and ¢/(1 + k). 


The terms describing the departure from equilibrium may now be calculated 
by retaining only the largest terms in all equations. Equations (20) and (21) become 


Ke, = c/(1 + k) (31) 
and equations (24), (22) and (26) consequently reduce to 
0 c 
1 08 r*\ 0 
0 dc @ 
By subtracting equation (34) from equation (33), we get 
1 0 0/ ke r2\ 0 
The solution of this equation, under condition (18), is 
1 d/ ke ] 
=: — ——. — 2s? + 6r? — 3— (3% 
— si + ‘) ( 8) 
The value of 5, as calculated from equations (32) and (36) is 
0 c s? 0/ ke 0 


The equation for c may in principle be found by substitution of equations (36) 
and (37) in equation (28). These calculations may, however, be greatly simplified 
by firstly estimating the relative importance of the various terms and discarding 
those that are negligible. We define a parameter by 

a = maximum [s?v/D, D/v, (Tv)/(1 + k)] 
Let f constitute a measure of the band width. c and its various derivatives will change 
from the maximum to the minimum value in a distance of the order of magnitude of f. 
This may be stated as 


p+aq p+q p+aq 
dx 4x B 


where the symbol ~ is used to denote an order of magnitude. The band will move 
along the column at a rate not differing much, if at all, from V. It then follows that 


+ + + 
at\axrats) Ath” ~ axPat 
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By a repeated application of equations (39) and (40) one finds 


Let ¢ denote any one of the quantities v, V, D, T and k. Under variable conditions, 
the variation in ¢ along the column length at constant time, and during the experi- 
ment at constant position, will normally be of the order of magnitude of the value 
of ¢ at any point. It now follows by the method used to determine equation (41) that 


tad, 1 \P/V 

The relation between a, 8, and L is best observed by introducing 

i.e. 

Our previous assumptions may be conveniently generalised a as 

y <1 (47) 


and the various orders of mugnitede may be denoted as powers a. y. Equation (47) 
is verified in paragraph 6, and our previous assumptions may now be shown to be 
included in equation (47). Equation (12) may be written as 


<a 


i.e. 

From equations (36) and (87), 


so that equation (29) is satisfied. 


Equation (28) ~ be symbolically written as 


where the order of sical of each bracketted term in equation (28) has been 
separately denoted. If the first two orders of magnitude are retained, it is clear that 
the second term on the right hand side of equation (28) may be neglected. In each of 
the other terms on the right hand side, the differentiation may be carried out accord- 
ing to the product rule, after insertion of equations (36) and (37). In every case, it 
follows from equations (41) and (42) that the term of lowest order will contain deriva- 
tives of c only, and not of v, V, k, D or T. The following example serves to demonstrate 
this point: 
0 ov We Ve. Ve Ve 
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In evaluating the terms on the right hand side of equation (28), v, V, k, D and T may 
therefore be regarded as constants. Equations (28), (36) and (37) now become 


dc 0 D 08 r?\08 

ot ax! 1 + k dx? ox s? *( (53) 
1 k @ r*\dc 

5, = ——| (3s? — 6r?) +- V[ —2s2 + 6r?2 — 3— |_— 54 

T @ k dc 


Substitution of equations (55) and (54) in equation (53) yields the general differential 
equation describing band development: 
dc ye ev. ( D 4 s*v? (1 + 3k)\0%c 
dt ox ax 1+k 48D 


ox 2 


(3) Solution of the general equation in terms of distance- time coordinates. As the 
total mass of solute will not affect the discussion, we impose the condition 


+4 00 
— 
on c. Let 
+00 
— 


designate the centre of mass of the solute band. In order to take into account the 
variation of v, V, D, k and T along the band, these quantities may be expanded 
‘in a series 


n=0 n! Ox"/z 


about z. The subscript z denotes the value of a quantity at the centre of mass, and it 
is important to note that it transforms all quantities to functions of time only. For the 
purpose of evaluating orders of magnitude, (x —z) may be replaced by f. By substi- 
tuting equation (59) in equation (56) and retaining terms of the first two orders of 
magnitude, we get 


ac aV\ Jac. /av 


(1+k)? 24D (1 +k)? 


c will vanish sufficiently strongly at large distances for the relation 
ax? 
to hold for all n and p. The above relation is verified by equation (76). On applying 


n 
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4 


the operator | xdx to all terms of equation (60), we find by partial integration, 
D 


and by using equation (61), 


dt i+k/. 


When v and k do not vary along the column length, k may be interpreted as the ratio 
of the average time spent by each molecule in the immobile phase to the time spent 
in the mobile phase, during the passage of the molecule through the capillary. Equa- 
tion (62) then shows that the average velocity of the solute molecule, while in the 
mobile phase, is equal to the average carrier gas velocity. Under variable conditions, 
the average transport velocity of the solute is seen to be given by the same expression, 
provided that the values of v and k at the centre of mass are used. 


Equation (60) may now be transformed to the centre of mass coordinates, defined 
by 


We shall use the symbols = and ~ to denote partial differentiation to time at con- 
1 


stant x and x,, respectively. The application of the operator = to c or its derivatives 
1 
yields a term of higher order of magnitude than that resulting by the application of 


= We may thus assume that 
dc 
at, 

Equation (64) is verified in paragraph 6. By the same arguments as before, this rela- 


tion will also apply to the various derivatives of c. By transforming equation (60) by 
means of equation (63), and retaining only terms of the lowest order, we find 


(64) 


2 
ot, Ox,/20x, Ox,/2 2 
where 
2 2 
2D 4 6k + 11k*) vs? (6) 
v 24(1+k)? D  (1+k)? 
Let c’ denote the r moment of c with respect to z, i.e. 
+00 
By partial integration, using equation (61), one finds 
+00 a+q 1 d@cp—a 
dtidx? (p — q)! dt} 
for p2>q 


for p<q 


ha 
D 
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+00 
On applying the operator | x? dx, to each term of equation (65), the result is 
— 


2 
dt Ox,/z 
when it is kept in mind that, from equation (57), 


and that the variance with respect to distance from the centre of mass is given as a 
function of time by 


ox 


Using equation (62), equation (69) may be written as 


dz Ox,/2 
The solution of equation (72), corresponding to the initial condition 

is given by 


° 2’\Ox 


Equation (65) has to be solved under the condition 


where 8(x,) denotes the Dirac delta function. Assuming an approximately Gaussian 
function, a solution in the form of an infinite series may in principle be obtained by 
evaluating the successive moments of c with respect to z.1* The Gaussian function 

2)—1/2 (x — z)? 
20% 
is, however, found to satisfy equation (65) identically when use is made of equation 
(69). 


(4) The form of the elution curve. The situation at the column exit is described by 


20? 

As the band width is, from equations (45) and (47), much smaller than the column 
length, o2 and V, may respectively be replaced by o2, and V, during the time interval 
in which the solute concentration at the column exit is significant. The subscript 0 
refers to the value of the various quantities at the centre of mass when the latter 
passes through the column exit. The relation between t and z is given, from equation 
(62), by 


4 
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The retention time is accordingly given by 


By subtracting equation (79) from equation (78) we get 
La,’ 
2V, Vo 


as a lowest order approximation. Substitution of equation (80) in equation (77) 
yields 


1 2\—1/2 (t — 
207 
where 


From equations (15), (30) and (81) the average concentration of solute in the 
gaseous phase at the column exit is given, in the general case where equation (57) 
is not applicable, by 


where 


ms*(1+k,)V, 
is a normalisation constant dependent upon the mass of the solute introduced into 
the column. 


By using equations (62), (74), (79), (82) and (83) a complete description may be 
given of the efficiency with which various separations may be carried out. When the 
values of the different parameters, e.g. v (x, t) and k (x, t), are known as functions 
of x and t, z(t), (the position of the centre of mass of the solute band at the instant t), 
and t(z), (the inverse function, i.e. the time taken by the centre of mass of the solute 
band to reach the point z), may be calculated by solving equation (62). This is best 
seen by writing equation (62) as 


dz z,t 
(85) 
dt 1+k(z,t) 
The retention time is then given by 


and the value of each quantity at the centre of mass may be found as a function of 
the position of the centre of mass, e.g. 


1 + k{z, t(z)] 
Once V, and h, are known, o; may be calculated from equations (74) and (82). The 
height equivalent per theoretical oni is then given by 


4 
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When the various operating conditions are functions of time only, as will approxi- 
mately be the case when the temperature of a relatively short, wide, capillary as a 
whole is varied during elution, it follows from equation (85) that 


t 

ol + k(t) 

from which the retention time may be determined. Using equations (74) and (82), 

is given by 


2 1 
= h,dz ee ee oe oe ee ee oe 90 
AR ( 


For techniques in which the operating factors are time-independent functions 
of x, t, and o7 may be directly obtained by integration of equations (74) and (79), 
using equation (82): 


L 
o 
L 2 
v*(x) 
In the derivation of equation (92) use has been made of the relation 
ox z dz 


which is valid for the case in which V is time independent. 


When a capillary is considered in which no variation of operating conditions 
occurs (this situation has been discussed by Golay*), the physical significance of 
h(x,t) becomes clear, as it follows from equations (88), (91) and (92) that h is then 
equal to H. The quantity h(x,t) of a column may thus be called the “local height 
equivalent per theoretical plate”, provided that the latter term is defined as the 
height equivalent per theoretical plate of a second, hypothetical, column in which 
the column conditions are constant and equal to conditions of the first column at the 
point x and the instant t. 


(5) The effect of the compressibility of the carrier gas at constant temperature. The 
pressure drop along the column is related at constant temperature to the average 
carrier gas velocity by 

s2 
and 

where the subscript i denotes the value of the various quantities at the column inlet. 


Since the diffusivity of the solute in the carrier gas may be taken?‘ to be in- 
versely proportional to the gas pressure, it follows that 
DP = D/P; oe oe ee oe oe oe (96) 


1960, Vol. XIII JOURNAL OF THE SOUTH AFRICAN CHEMICAL INSTITUTE 109 


Introducing the ratio of inlet to outlet pressure, 
it follows from equations (94) and (95) that 


The local plate height given by equation (66) may, from equations (95) and (96), 
be written as 

where A and B are constant pry given aby 


24(1 +k)? D, 
For constant k, we find from oqentions (88), (91) and (92), 


Substitution of equations (98) and (99) into equation (102) yields, after the integration 
is performed, 


(100) 


H = Ay,(p) + Buy,(p) 
where 
_ 9(p* — 1)(p* — 1) 
_ 3 p(p? — 1) 


These equations are in agreement with the result derived by Giddings et al.5.® and 
may be compared to the equation derived by Littlewood,’ 


H = A + Buy, (p) 
where 

and the equation given by De Wet, ® 

where 

ys(p) = (8Sp—3)/(5p) .. .. .. . a6 


The functions of p given by equations (104), (108), (107) a and (109) have been plotted 
against 1/p in Fig. 1. Equations (106) and (108) do not contain any correction for A. 
The correction for B is given within 2% by equation (108) for p<5, with a maximum 
error of 6-6% for p= ©, and is given by equation (106) within 6% only for p<2-5, 
with a maximum error of 333% for p= ©. The total maximum error possible 
is 111% with equation (108) and 33-3% with equation (106). Equation (108) appears 
therefore to be a more accurate approximation than equation (106). It may further 
be noted that the height equivalent per theoretical plate will fluctuate between 1 
and 1-5 times the value of the local plate height at the column inlet. 


| 
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° 
N 
ys(p) 
a 
y2(p) 
| 
yi(p) 
0 4 8 10 


Ip 


Fig. 1.—y,(p), y2(p), y4(p) and y,(p), respectively calculated from equations 
(104), (105), (107) and (109), plotted against 1/p. 


As a measure of the efficiency of separation of a mixture of solutes, the height 
equivalent per theoretical plate does not take into account the influence of the reten- 
tion of the solutes. This point is illustrated by a column in which k = o, for which H 
will have a finite value, but in which the separation of a solute mixture is impossible. 
This defect may be corrected? by the introduction of a height equivalent per theor- 
etical plate based on the apparent retention time, 


1 + k)? 


k?2 


(110) 


€ 

| 
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which should be minimised for optimum separations. By introducing the ratio of 
the liquid film thickness to the ee radius 

and by using equations (1), (11), (100), 101), (103) and 1 (110), the full expression 
for Hgar in terms of the inlet conditions is 


_ (Di (1 + 6k + 
Hgar = (1 — 0)? (p) 
+ (112) 


From equations (95) and (96), Hgap may be expressed in terms of the outlet 
conditions: 


D, (1 k)? 6k + 11k? 
2 
where 
3 (p* — 1) 


y3(p) is similar to the factor derived by James and Martin?® to correct for the effect 
of the pressure drop on retention volume, and has been tabulated!’ for 1 < p < 3. 


Hgag may also be expressed in terms of the average carrier gas velocity, 
obtained from equation (98) as 


L 


It is of interest to calculate the ultimate efficiency theoretically possible in a 
coated capillary column. This will theoretically be realised when D, = © and y = 0. 
Equation (112) may then be minimised with respect to v,, to yield 


(1 + = + 6k + =) 
H =2 ‘ —@ 
( BAR) min k 12k2 ) ( ) 
This expression will be a minimum for k = ™, 

(Hpar)min = 2(11/12)*/#r, (1 — 8) ee 


The factor (1 — @) accounts for the reduction of the radius of the mobile phase by 
the stationary phase. The minimum H is given under the same conditions by 


2 1/2 
12(1 + k)? 
with an absolute minimum value, for k = 0, of 
(H) min = 2(1/12) 2/ 2r,(1 — 0) os 


At this point separation will, however, not be possible, which further illustrates the 
unsuitability of H as an efficiency factor. 


§ 
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(6) Discussion of assumptions. The assumptions made in equations (47) and (64) 
may now be verified. From equations (38) and (66) it is seen that 


Differentiation of equation (76) with respect to time at constant x, yields 

dc 1 \de? 

ot; 
From equations (120) and (69), 

dt 


As distances x, much larger than £ are not of any importance, equations (121) and 
(122) reduce to 


aV V 
at, 
which thus verifies equation (64). 


In order to verify equation (47), 8 may be chosen as equal to o,. Equation (74) 
then becomes 


B®? ~hz ee (124) 
From equations (43), (120) and (124), 

y ~ maximum [(h/z) 1/2, ~ (h/z) 1/2 
Equation (47) then becomes 


When the pressure ratio p is not very large, it follows from equations (98), (99) and 
(103) that 


which reduces equation (126) to 


where N is the total number of equivalent theoretical plates. For a column with a 
large plate number, the approximation made in equation (47) will thus be valid for 
all regions except a relatively short section near the column inlet, the length of which 
is of the order of magnitude of H, and the influence of which will become zero as N 
tends to infinity. When p is large, it follows from equation (98) and (99) that equation 
(127), as a result of the resistance of the solute to mass transfer in the stationary 
phase, is not applicable in the exit region, and should be replaced in the extreme 
case where p = ©, by 


The difference between equations (127) and (129) becomes marked when 

It then follows from equations (129) and (130) that 


(130) 
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from which it may be seen that deviations from equation (127) will occur, for large N, 
in a relatively short section near the column exit. As previously pointed out, this 
effect is due to resistance of the solute to mass transfer in the stationary phase. The 
maximum deviation 4o, caused by the column exit may approximated from equation 
(10) by 

From equations (88), (99), (101) and (115) it follows that 


The influence of the region in which equation (127) is not applicable is thus negligible 
for large N. 


(132) 


On the basis of the foregoing discussion it may be concluded that the assumptions 
which have been made in the derivation of the general equations will not affect the 
accuracy of the equations provided only that the number of theoretical plates is large. 
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LIST OF SYMBOLS 


Symbol Dimensions 

A ] .. Convenient factor defined by equation (100) 

B t .. Convenient factor defined by equation (101) 

C ml~“*t .. Normalisation constant dependent on total mass of 
solute 

c(x,t) ml-* .. Convenient factor defined by equation (13) 

c(t) ml"? .. 1" moment of c with respect to position 

c,(x,r,t) ml .. Concentration of the solute in the mobile phase 

C(x, t) Radial average of c, 

C,(x,T,t) ml .. Concentration of the solute in the stationary phase 

.. Radial average of c, 

D(x,t) 1% .. Diffusivity of the solute in the mobile phase 

D, (x,t) .. Diffusivity of the solute in the immobile phase 


H ] .. Height equivalent per theoretical plate 


2 
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Height equivalent per theoretical plate based on 
apparent retention 

Local height equivalent per theoretical plate 

Subscript designating values at the column inlet, at 
constant temperature 

Mass flux of the solute in the mobile phase 

Mass flux of the solute from the mobile to the immobile 
phase 

Ratio of the solute concentration in the mobile phase 
to the concentration in the immobile phase, at 
equilibrium 

Ratio of the mass of solute in the immobile phase to 
the mass in the mobile phase, at equilibrium 

Column length 

Total mass of solute 

Number of equivalent theoretical plates 


Subscript denoting values at the centre of mass of the 
solute when the centre of mass passes through the 
column outlet 

Gas pressure at constant temperature 

Ratio of the inlet to the outlet pressure at constant 
temperature 


Ratio of the volume of the stationary phase to the 
volume of the mobile phase 


Radial distance from the centre of the column 

Radius of capillary 

Radius of mobile section 

Convenient factor defined by equation (11) 

Time as measured from instant of injection of solute 

Retention time 

Carrier gas velocity 

Velocity of the centre of mass of the solute when the 
latter is at (x,t) ; 

Radial average of the carrier gas velocity 


Longitudinal average of the carrier gas velocity, as 
defined in equation (115) 


Thickness of the retentive layer 
Distance from column inlet 


Hgar l 
h(x,t) 
i None 
j(x.r,t) 
js(x,t) 
K(x,t) None 
k(x,t) None 
M m 
N None 
None 
P(x) 
Pp None 
s 
T(x,t) t 
t,t, t 
t, t 
u(x,r,t) 
V(x,t) 
u(x,t) 
x a 


1960, Vol. XIII JOURNAL OF THE SOUTH AFRICAN CHEMICAL INSTITUTE 115 


Symbol Dimensions 

Xy l .. Distance from centre of mass of sample 

y(p) None .. Various functions of the pressure ratio, defined in 
equations (104), (105), (107), (109) and (114) 

z(t) l .. Position of the centre of mass of the sample. The 
subscript z denotes the value of a quantity at the 
centre of mass 

a(x,t) l .. Convenient factor defined by equation (38) 

B(x,t) 1 .. Order of magnitude of the band width 

y(x,t) None .. Convenient factor defined by equation (43) 

5,(x,r,t) ml? .. Convenient factor defined by equation (16) 

8,(x,t) ml? .. Convenient factor defined by equation (17) 

n ml~“'t—' .. Carrier gas viscosity at constant temperature 

6 None .. Ratio of the thickness of the retentive layer to the 
radius of the capillary 

p(x,t) ml! .. Total mass of solute per unit length of the column 

oO, t .. Standard deviation of the concentration of the solute 
with respect to time at the column outlet 

o,,(t) l .. Standard deviation of the concentration of the solute 

with respect to distance 

$(x,t) Convenient factor designating v, V, D, T or k 
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GAS-LIQUID CHROMATOGRAPHY IN CAPILLARY COLUMNS. PART 2. 

FACTORS AFFECTING THE COLUMN EFFICIENCY: THE DISTRIBU- 

TION COEFFICIENT, THE DIFFUSIVITY OF THE SOLUTE IN THE 

GASEOUS AND LIQUID PHASES, AND THE THICKNESS OF THE 
LIQUID LAYER 


by 
P. C. HAARHOFF and VICTOR PRETORIUS 


OPSOMMING 


Die invloed van die distribusiekoéffisiént, die diffusiwiteit van die opgeloste stof in die gas- 
en vloeistof-fases, en die dikte van die vloeistoflaag op die effektiwiteit van ’n kapillére kolom is 
teoreties ondersoek. Die hoogte ekwiwalent aan ’n teoretiese plaat, gebaseer op skynbare retensie, 
word gebruik as ’n maat van die effektiwiteit. Grafiese resultate word gegee waaruit ’n skatting 
gemaak kan word van die effektiwiteit van die kolom wanneer die waardes van die parameters 
hierbo genoem, bekend is. Die bestaan van ’n optimumwaarde van die dikte van die vloeistoflaag 
word aangetoon. Dit blyk dat baie hoé kolomeffektiwiteite moontlik is wanneer beide die dis- 
tribusiekoéffisiént en die verhouding van die diffusiwiteit van die opgeloste stof in die gasfase tot 
dié in die vloeistof-fase relatief klein is. 


SUMMARY 


The effect of the distribution coefficient, the diffusivity of the solute in the gaseous and liquid 
phases, and the thickness of the liquid layer on the efficiency of a capillary column has been 
theoretically examined, using the height equivalent per theoretical plate based on apparent 
retention as a measure of the efficiency. Graphical results are given which may be used to estimate 
the column efficiency when the values of the parameters listed above are known. The existence 
of an optimum value of the thickness of the liquid layer is demonstrated. It appears that very high 
column efficiencies are possible when both the distribution coefficient and the ratio of the diffusi- 
vity of the solute in the gaseous phase to that in the liquid phase, are relatively small. 


The efficiency with which chemical mixtures may be separated by gas-liquid 
chromatography in coated capillary columns is dependent on a number of factors. 
The following expression relating these factors to the column efficiency has been 
derived?: 


(1 + k)?D, 
H = 1—6)2_: 

BAR [ k? 24k? ) ) D, (P) 

2 6? r? 

+ Ye (P) (1) 
where 

D, = diffusivity of the solute in the gaseous phase at the column inlet 
D, = diffusivity of the solute in the liquid phase 


Hgar = height equivalent per theoretical plate based on apparent retention ? 
K == ratio of the concentration of the solute in the gaseous phase to the 


concentration in the liquid phase, at equilibrium 


St 


1960, Vol. XIII JOURNAL OF THE SOUTH AFRICAN CHEMICAL INSTITUTE 117 


k = ratio of the mass of the solute in the liquid phase to the mass in the 


gaseous phase, at equilibrium. 
P = ratio of the inlet to the outlet pressure 
5 I, = radius of the capillary 
, Vi = radial average of the carrier gas velocity at the column inlet 
6 = ratio of the thickness of the liquid layer to the capillary radius 
9 (p* — 1) (p* — 1) 
(P) 8 — 1)? ( 
3 p(p* — 1) 
Ya (P) = 5 1) (3) 
20 — 6? 
(4) 
(l1— 6)? K 


Equation (1) may be used to study the way in which the column efficiency will 
change under a variety of experimental conditions and also to determine the experi- 
mental conditions which should be used to give the best possible separation of a 
particular mixture. The best method of doing this is to assign to each parameter in 
equation (1) a value chosen from the range of possible values which might be en- 
countered in practice, and to calculate the column efficiency in each case. It is obvious 
from equation (1) that a very large number of sets of values of the parameters is 
is possible. The method is further complicated by the fact that many of the parameters 
e are mutually interdependent so that the value of one parameter cannot, in general, 
be changed without affecting the values of one or more of the other parameters at 
the same time. The relationship between mutually dependent parameters is also 
frequently complex and sometimes even unknown. A complete description of the 
effect, on the column of efficiency, of every possible value of each parameter, would 
obviously become excessively laborious. 


In order to reduce the number of calculations, whilst at the same time retaining 
the widest possible range of parametric values, the parameters may be treated as 
dimensionless groups, as follows: 


We define 
Y = Hgar/ro oe oe ee ee ee ee (5) 
X = v,1,/D, (6) 
4 Substitution of equations (5) and (6) in equation (1) yields 
Y = (A/X + BX) y, (p) + CXyz (p) 
where 
2 _ 2 
24k? 
2 
(10) 
3 kD, 
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Equation (7) is simpler for the present purpose than is equation (1) in that values 
only of X and not of v;, r, and D, separately, need be considered. A further simplifica- 
tion of equation (7) may be made by ignoring the pressure correction terms, viz. 


It follows from equations (2) and (3) that the first two terms in equation (11) can 
never differ by more than 11-1% from the corresponding terms in equation (7), and 
the third term never by more than 33-3°%. Equation (11) thus gives a precise descrip- 
tion of the processes in short, relatively wide-bore capillary columns as well as provid- 
ing a description of the general trends in long, relatively narrow capillary columns. 
Y is given by equation (11) as a function of only four independent variables, i.e. 
X, K, 6 and D,/D,. A comprehensive study of the effect of various values of these 
parameters on Hxg,p is now possible. It should be noted that this is only true when 
the equation for Hg,p is expressed in terms of the values of the parameters which 
pertain to the column inlet. 


It follows from equations (5), (6) and (11) that, at constant K, @, D,;, D, and r,, 
there is an optimum value of the carrier gas velocity 


(X)opt = (Vi = WVA/(B + C) 
for which Hgap will be a minimum, 


As D, is dependent on the inlet column pressure, the optimum value of v; is obtained 
by varying the gas pressure at the column outlet. From equation (13) it may be seen 
that the minimum value of Hg,, is directly proportional to r, when K, @ and D,/D, 
are kept constant, and is independent of the absolute values of the diffusivities of 
the solute in the mobile and stationary phases. Higher efficiencies will thus be realised 
with narrow than with wide capillary columns. The optimum flow velocity of the 
carrier gas is proportional to 1/r, when K, @, D; and D, are kept constant, and narrow 
capillaries are thus further desirable from the point of view of rapid analyses. It 
follows from the above discussion that the effects on Hap of varying K, 6, D, and D, 
at constant r), and of varying r, at constant K, 6, D; and D,, may be superimposed 
upon each other, and r, may thus be regarded as a constant in the work which follows. 

The following hypothetical values of the parameters are representative of the 
range of possible values likely to occur in practice: 


K = 0-001, 0-01, 0-1, 1 
D,/D, = 103, 105,5 x 105 
6 =0-002 to 0-3. 


Equations (8), (9), (10) and (11), and values of the parameters chosen from the above 
list, have been used to calculate Y for various values of X. The results are shown in 
Figs. 1-3. The minimum value of Y for each value of K and D,/D, listed above has 
been calculated from equation (13) for various values of @ (Fig. 4). 


The results in Figs. 1-3 show that the optimum carrier gas velocity increases 
as the resistance of the solute to mass transfer in the liquid phase is lowered, at 


: 
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| 

4 
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constant K and D;, by decreasing the percentage of liquid and/or by increasing the 
diffusivity of the solute in the liquid phase. This effect is marked in all cases except 
when K and D,/D, are both very small and may be explained as follows. It is apparent 
from equations (11) and (12) that, when the ratio of (B + C) to A is small, the in- 
fluence of the solute resistance to mass transfer (BX and CX) becomes of the same 
order of magnitude as that of the longitudinal diffusion (A/X) only at relatively high 
carrier gas velocities. Equations (10) and (12) show that high optimum carrier gas 
velocities also result for large values of the diffusivity D; of the solute in the mobile 
phase; the minimum value of Hy, p, however, becomes larger when D,/D, is increased 
(Fig. 4). The time taken for an analysis may thus be more profitably reduced by 
employing large values of D, than by increasing D;. When K is decreased while @, 
D; and D, are kept constant, the optimum carrier gas velocity decreases slowly for 
small values of @ and D,/D,, while the opposite is true for larger values of @ and D,/D,. 


The time for analysis may also be shortened by increasing the gas velocity 
beyond the optimum value. However, it is evident from Figs. 1-3 that the column 
efficiency is at the same time decreased. This decrease may be determined by expres- 
sing Hg,r and vy, in terms of their optimum values, as follows: Let 


and v, = Vi/(Vi)opt wi (15) 
= 
= FIG. 1—The influence of the inlet gas 
velocity on the height equivalent per 
theoretical plate based on apparent 
retention, for K = 0-001 and the 
an following values of @ and Dj/D,: 
a 6 
~ 0-002 0-01 0-05 
« 
3 
3 5x105| a, | b | c | 
Di/D, 10° | a b, | 
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\ 10? | a, | & 
ak 
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a FIG. 2—The influence of the inlet gas 


} velocity on the height equivalent per 
theoretical plate based on apparent 
retention, for K = 0-01 and the fol- 


2 lowing values of @ and Dj/D,: 
a 6 
0-002 0-01 0-05 
| 5x a, | by | | 
Di/D, 10° | a, | be | Cc, | 
103 | a, b, | Cy | 


Y= Hean/re 


1 afi 1 


. 2 4 6 8 10 


X= vi re/Di 


From equations (11), (12), (13), (14) and (15) a particularly simple equation is found, 
viz.: 

(Hpar)q = (1/2) (I/vq + Va) 
from which it ‘ion that an increase of the carrier gas wdlolty by a certain fraction 
of the optimum carrier gas velocity will necessarily result in an increase of Hgapr by. 
a definite fraction of (Hgar)opt y- Once a high optimum carrier gas velocity is found 
with a satisfactory low value of Hgar, it may thus be increased by a correspondingly 
large value for a stipulated percentage of loss in efficiency. The precise relationship 
follows from equation (16). It may be noted that (Hg,r), changes slowly from 1 to 
1-25 when v, changes from 1 to 2, while a more rapid rate of increase 
(d(Hgar)q/dv, = 0°5) is found for high values of v,. 


According to Fig. 4, Hg,pr is, in addition to r,, dependent on 6, K, and D,/D,, 
but not on the absolute values of D; and D,. If (Y),,.x is plotted against @, three 
distinct regions occur. For k < 2, i.e. @ S K, both A and B increase rapidly as the 
liquid film thickness is made smaller, since the solutes are then not sufficiently 
retarded by the liquid for good resolution of the solute bands. From equations (4) 
and 3 C may be approximated by 

= (1/3) 6K (D,/D,) 
which it may be seen that C ter: the under discussion. 
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This is illustrated in Fig. 4, as (Y),,¢ x is relatively insensitive to variations in D,/D, 
for small values of 0, especially for K = 0-001. For 6 2 3K and D,/D, & 104, the 
resistance of the solute to mass transfer is located almost entirely in the liquid phase, 
and increases when more liquid is used. An intermediate region, K < 6 < 3K, may 
be distinguished, in which (Y).,,. x varies very slowly with @, and in which most of 
the curves in Fig. 4 exhibit an optimum percentage of liquid. The column efficiency 
can thus not be decreased beyond a certain limit by decreasing the thickness of the 
liquid layer, and consequently the resistance of the solute to mass transfer in the 
liquid phase. This result is somewhat surprising, since the equation derived by Golay *, 
in which the efficiency is expressed in terms of the height equivalent per theoretical 
plate, predicts that the column efficiency increases continuously as the thickness of 
the liquid layer is decreased, once k < 1. The reason for this discrepancy lies in the 
employment of the height equivalent per theoretical plate as a measure of the column 
efficiency, instead of the more correct height equivalent per theoretical plate based on 
apparent retention. In Fig. 5 the minimum value (Y),,. x,,, Of (Y)opt x in the region 
0-002 < @ < 0-3 has been plotted against D,/D, for K = 0-001, 0-01, 0-1 and 1. 


Very high column efficiencies may be obtained from small values of K and D,/D,. 
The smallest Hg, p-value possible has been shown to be? 


240 


FIG. 3—The influence of the inlet gas 
velocity on the height equivalent per 
theoretical plate based on apparent 
retention, for K = 0-1 and the fol- 
lowing values of @ and Dj/D,: 
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7 FIG. 4—The influence of the thickness 

ss of the liquid layer on the minimum 
height equivalent per theoretical plate 
based on apparent retention, for the 
following values of K and Dj/D,: 

K 

8 0-001 0-01 0-01 1 
5 x 105 | bs | cs | as | 
Di/D, 10% | ag | | | dy | 

oe} a | | | a, | 
2 
0-002 0-01 0-3 


for K = 0, D,/D, = 0. The factor ( 1 — @) accounts for the reduction of the radius of 
the mobile phase by the stationary phase, but it would be impractical to reduce 
Hyar by this means. The ultimate efficiency possible may consequently be taken as 

(Hgar/To)min = 1°915 (19) 
Extensive experimental data on the eficiency of capillary ooluene have not, as yet, 
appeared in the literature. An analysis of the results given by Golay® for the elution 
of o-xylene with helium as carrier gas and didecyl phthalate as stationary phase, 
shows that with 6 = 0-005 and K = 0-00074, a value Hgar/r, = 83 was obtained, 
which approaches the best possible value within a factor 4-3. 


It is clear from Figs. 4 and 5 that K and D,/D, should be made as small as possible 
in order to facilitate favourable performance. This may best be achieved by using a 
liquid which interacts strongly with the solute, and has a low viscosity‘ in order to 
increase the diffusivity of the solute in the stationary phase. The use of elevated 
pressures and dense carrier gases will decrease D,,5 but, as previously pointed out, 
will result in slow optimum flow speeds. For D,/D, < 104 and K < 0-01 or for K < 
0-001 and a wide range of values of D,/D,, it becomes unnecessary to reduce K or 
D,/D, any further, as this would not result in a large gain in efficiéncy. In this region, 
very thin liquid layers should be used to facilitate rapid analyses. When weakly 
interacting substances are analysed, i.e. K 2 0-1, useful results will be obtained only 
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for K ~ 0-1 and D,/D, < 104, provided that a relatively large percentage of liquid 
is used. As a result of insufficient retention, solutes for which K ~ 1 yield impossibly 
high values for (Hg,r/T,)opt v» even When values of @ as large as 0-3 are used. 


In determining the best value of 6 for a given value of K and D,/D,, column 
efficiency as well as speed of analysis should be taken into account. The values of 8 
indicated by the circles in Fig. 4 have been chosen from these considerations. The 
value of the velocity of the solute band at the minimum value of Hgar, (Vi)ope/(1 + k), 
multiplied by r,/D, in order to obtain a dimensionless factor, has been calculated for 
each point from equations (4) and (12). In Fig. 6 a plot is given of (v;),»:To/D; (1 + k) 
against K for various values of D,/D,. 


It is seen that for K ~ 1, the optimum band velocity is very low, and Hgap is 
very large. Under such conditions, there is little point in striving to obtain a reason- 
able efficiency: although this could be done by using extremely long columns and thick 
liquid layers, the time required to elute the solute would become excessive. For 
reasonably large values of D,/D,, strongly interactive solutes are eluted more slowly 
than are solutes which interact moderately with the solvent. This tendency would 


(Y opt = (Haar/to)opt 


1o* 10° 
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FIG. 5—The influence of Di/D, on the minimum height equivalent per theoretical 


plate based on apparent retention, at the optimum value of @ in Fig. 4, for the 
following values of K: (a); K = 0-001, (b); K = 0-01, (c); K = 0-1, (d); K = 1 
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0-001 0-01 0-1 


K 


FIG. 6—The influence of the distribution coefficient on the value of the optimum 
band velocity, (vi)opt/(1 + k), at the points indicated on Fig. 4, for the following 
values of Dj/D,: (a); Di/D, = 10%, (b); Di/D, = 105, (c); Di/D, = 5 x 105 


08 12 16 20 


CVidopt re/Dt (1+k) 


0-4 


be offset by employing values of @ smaller than the value 0-002, which has been chosen 
as the lower limit of @ in the present study. As a lower limit of @ will exist in practice, 
it would appear that the use of very low K values, when D,/D, 2 105, would have 
the disadvantage of increasing the time of analysis without at the same time yielding 
a significant improvement in the column efficiency. A more satisfactory balance 
between these factors should thus then be obtained by using K values of roughly 
0-003. 


One of us (P. C. H.) is indebted to the South African Atomic Energy Board for 
the award of a bursary, and the other (V. P.) to Messrs. African Explosives and 
Chemical Industries, Ltd., and to the Council for Scientific and Industrial Research 
for research grants. ; 


Department of Physical Chemistry, 
University of Pretoria. Received August 16, 1960 


REFERENCES 


1 Haarhoff and Pretorius, preceding paper. 

2? Haarhoff and Pretorius, to be published. 

3 Golay, Gas Chromatography, edited by Desty, Butterworths Scientific Publications, London, 
1958, p. 36. 

# Othmer and Thakar, Ind. Eng. Chem., 1953, 45, 589. 

3 Gilliland, ibid., 1934, 26, 681. 


é 

a 

7 

4 


A 


1960, Vol. XIII JOURNAL OF THE SOUTH AFRICAN CHEMICAL INSTITUTE 125 


NOTES 


THE USE OF SMALL QUANTITIES OF SULPHURIC ACID AND 
POLYPHOSPHORIC ACID IN THE VON PECHMANN REACTION 


by 


B. E. BARRIS and S. S. ISRAELSTAM 


Greig Gass ' showed that the relatively large quantities of concentrated sulphuric 
acid used by von Pechmann? in his synthesis of hydroxycoumarins by the interaction 
of phenols and f-keto esters were unnecessary, and that only a few drops of the acid 
were sufficient, provided that acetic anhydride was added to take up the water 
formed in the reaction. 


Some years ago it was found that the addition of the acetic anhydride was 
unnecessary and that a relatively small quantity of sulphuric acid alone was quite 
sufficient to give good yields of 7-hydroxy-4-methylcoumarin by heating together 
fifth molar quantities of resorcinol, ethyl acetoacetate and 0.1 g concentrated sul- 
phuric acid for about 30-40 minutes at 120°C. Recently the opportunity presented 
itself to continue with this investigation and the formation of hydroxycoumarins 
from various phenols and f-keto esters using small quantities of sulphuric acid has 
been studied in detail. Furthermore, the use of polyphosphoric acid was also investi- 
gated. Koo has shown that this reagent has certain advantages over sulphuric acid 
and can be used as a condensing agent in this reaction, but as with von Pechmann’s 
use of large quantities of sulphuric acid, he used large quantities of polyphosphoric 
acid. Thus, for tenth molar quantities of resorcinol and ethyl acetoacetate as much 
as 160 g of polyphosphoric acid were used. We have found that comparable yields of 
coumarins are obtained using very much smaller quantities of polyphosphoric acid 
but longer heating is necessary than in the case of sulphuric acid. In general, better 
yields were obtained with sulphuric acid. 


In order to find the optimum conditions for this reaction, a series of experiments 
was carried out, varying the quantity of condensing agent and the time of the reaction 
in each case. These results are summarised in the Table. 


It will be noted that very good yields of coumarin were obtained in the case 
where “reactive” phenols (e.g. resorcinol) and “‘reactive’’ B-keto esters (e.g. ethyl 
acetoacetate) were condensed together, but that low yields of coumarin resulted 
when less “‘reactive’’ phenols (e.g. orcinol) and f-keto esters (e.g. ethyl a-methyl- 
acetoacetate) were used. In most cases the optimum temperature was found to be 
120-130°C and the time of heating varied from 5 minutes to 30 minutes, depending 
on the weight and reactivity of the phenol and f-keto ester. 


After this work had been completed, our attention was drawn to a patent by 
Boese* who had developed a method for the preparation of hydroxycoumarins also 
using small quantities of concentrated sulphuric acid, but using toluene as a solvent 
and refluxing the mixture for an hour in order to remove the alcohol and water 
formed in the reaction as an azeotropic mixture. Toluene is not required on the 
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laboratory scale. Furthermore, Boese did not list polyphosphoric acid among the 
condensing agents given in his patent. 

The method described here has been used for some years by students for the 
laboratory preparation of 7-hydroxy-4-methylcoumarin and a preliminary note has 
already appeared ®. 

Consideration has been given to the mechanism of the von Pechmann reaction 
in the light of these and other experiments and the results will be published elsewhere. 


EXPERIMENTAL 
Hydroxycoumarins. The following was the general procedure adopted: 
Equimolecular quantities of the phenol and the f-keto ester were heated 

together with various quantities of condensing agent in an oil bath at 120°-130°C. 

When the reaction appeared to be complete, i.e., when either the product solidified 

or when no further effervescence occurred, water was added, the mixture cooled, 

the crude product filtered off, washed with hot water and then with a small volume 
of 5 per cent. ethanol. The coumarin was crystallised from 50 per cent. ethanol. 


Where the product did not solidify ‘on the addition of water, the following 
methods were adopted: 

(a) If the product became pasty on the addition of water as in the case of 5,7- 
dihydroxy-4-methyl-3-ethylcoumarin, the coumarin was extracted by boiling the 
paste with water and allowing the solution to crystallise. This process was repeated 
a number of times. 

(b) If an oil was formed on the addition of water and did not solidify on cooling 
ethanol was added drop by drop until solidification occurred. 


The results of a number of experiments are summarised in the Table. Mixed 
melting points of the various coumarins with authentic specimens in each case 
showed no depression. The melting points of the crude coumarins in all cases were 
2°-6°C below that recorded for the pure material. 


We thank the Council for Scientific and Industrial Research for an award of an 
assistantship to one of us (S.S.I.). 
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VIR DIE AANDAG VAN OUTEURS 


Die Joernaal van die Suid-Afrikaanse Chemiese Instituut word gewy aan die publikasie van 
oorspronklike navorsing in alle vertakkings van chemie en chemiese teguologie, 

Alle bydraes moet in triplikaat ingehandig word; indien daar grafieke is moet een stel nou- 
keurig met Oosindiese ink afgetrek word op gladde wit Bristol-blad of tehenpapier, of op aftrekdoeh— 
die ander twee stelle kan afdrukke wees, of met potlood geteken word. 

Artikels moet so saaklik moontlik wees, verenigbaar met duidelikheid, en moet nie onnodige 
eksperimentele besonderhede bevat nie. Artikels moet in tékshkrif wees—dubbel spasiéring—net 
aan een kant van die papier, met ’n behoorlike kantruimte aan die linkerkant. Die papier moet 
nie-absorberend wees. 

Die buitemate van plate—insluitende die ruimte vir byskrif—moet hoogstens 14 duim by 10 
duim wees, en liefs nie minder as 7 duim by 5 duim nie. Die verhouding 1.4:1.0 moet in elk geval 
nagekom word. Byskrif en letters—hetsy in die kantruimte of in die tekening self—moet nie met 
ink gedoen word nie, maar liggies met gewone potlood of uitueebare bloupotlood, 

Elke artikel moet voorafgegaan word deur 'n opsomming, normaalweg van 50 tot 250 woorde 
lank, in elk van die offisiéle landstale. Die opsomming wat tweede verskyn, moet in dieselfde taal 
as die artikel wees. 

Verwysings moet in numeriese opeenvolging aan die end verskyn, in die volgorde; skrywer 
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Piping for the Highveld 
Power Station manufac- 


We are specialists in the 
design, manufacture, 
manipulation and. erec- 
tion of high pressure ther- 
mal and other pipework 
associated with Steam 
Power Stations and In- 
dustrial Process Plants. 


Expert planning plus the use of thoroughly reliable 
equipment will ensure your installation being both 
efficient and economical. This is the complete service 
we offer you. Moreover, we can place at your disposal 
not only the technical resources and equipment of 
our own organisation, but also those of world 
renowned overseas manufacturers with whom we are 
associated. We supply Stewarts and Lloyds Steel and 
Copper Tubes and Fittings: 

Hopkinson’s Steam Valves and Boiler Mountings: 
Glenfield & Kennedy Valves for water: 
Worthington-Simpson Pumps & Compressors: 


STEWARTS AND LLOYDS 


OF SOUTH AFRICA LIMITED 


YOUR PARTNER IN. PROGRESS 
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